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PREFACE 


The year of 1959 marked an epoch to be remembered by the coastal engineers 
in Japan for a new record of storm surge damage took place in the central part 
of the Japanese coast. The significance of coast protection was again manifested 
to this country where the coastal regions are highly utilized for industries and 
other productive activities as well as habitation. Thus numerous new subjects 
have been added to the field of coastal engineering which has only a short history 
in this country. 

Our intention in publishing the series of "Coastal Engineering in Japan" is 
to promote the researches and practices of this field of engineering through 
discussions from overseas and to reinforce exchange of knowledge between Japan 
and foreign countries. 


This volume contains the papers published in 1957 to 1958. 


Masashi Hom-ma 


Chairman 
Committee of Coastal Engineering 
November, 1959 Japan Society of Civil Engineers 


Tokyo, Japan 
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ON THE OCEAN WAVES DUE TO TROPICAL CYCLONES 


Sanae Unoki, Dr. Sc., Oceanographical Laboratory, 


Meteorological Research Institute, Tokyo 


1. Preface 

An attempt to clarify the structure of ocean waves in tropical cyclones such as typhoons 
and hurricanes has seldom been made because it is extremely difficult to take an exact and 
systematic wave-observation over an extensive sea covered by mountaineous waves. In the 
| meanwhile, the theory of ocean wave forecast has made considerable progress recently, but 
it is not yet satisfactorily powerful for ovr problems, in which the speed and direction of 
wind change rapidly with time and place. 

Under such circumstances, the present author has investigated the distribution of wave 
height, period and other elements in typhoons on the basis of plentiful wave materials 
obtained by visual observations from weather ships and military planes for about one hundred 
typhoons. In this paper, however, the outline of the results will be given since the full 
discussion was already made in other journals (see references). Moreover, it is hoped that 
these results will be amended with further accumulation of exact data, because both quality 
and quantity of our wave materials remain still insufficient, particularly in the vicinity 


| of the storm center. The wave height and period observed are those of prevailing waves 


which may correspond roughly to significant waves. 


2. Distribution of wave height 


Fig-l shows an example of mean distribution 


| of wave heights. The most remarkable feature is 
the asymmetry of height - distribution to the 
storm center, that is, the wave height is great- 
est in the right rear quadrant of the storm as 
pointed out by Cline (1926) and other investiga- 


| tors and smallest in the left front quadrant. 


This is attributed to the following conditions: 


| (1) Winds blow with greater velocity in the 
Fig-1 Mean distribution map of 
wave height (in meter) due to a 


2 typhoon whose central pressure 
or navigable semi-circle. (2) In the right semi- oi Co6G mb. 2 


right or dangerous semi-circle than in the left 


circle waves move in a direction near to that of 


movement of the storm, therefore, the fetch and 
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duration of winds become longer than in the left semi- 
circle. (3) In the rear semi-circle we can usually 
find remarkable swells which were originated in the 
front semi-circle and were left behind, consequently, 
wind waves generated in the rear are superposed on and 
raised by these swells. 

On the other hand, the distribution of waves 
depends remarkably on the speed of typhoon (C). Here- 


after, the term "wave velocity" is always used in a 


sense of "group velocity" because it could be ascer- 
tained that swells sent out from a typhoon travel with (Q) C:0~20 Klee 
the group velocity of deep water wave as theory demands. 
Jhen the storm is slower than the usual wave (Fig-2a, 
C<20 km/hr), the area of rough sea is wider in the 
right semi-circle than in the left semi-circle. On the 


contrary, when the storm is faster than the wave (Fig-2c, 


600SH 


C>40 km/hr), the wave is remarkably rougher in the 

rear semi-circle than in the front semi-circle since 

the wave arisen in the latter region is left in the for- 
mer region. When tke velocity of storm is near to that 


of wave (Fig-2b, C: 20~40 km/hr) the wave develops most 


conspicuously on account of the increase of duration anc 
fetch of wind, and also the asymmetry of wave-height 
distribution diminishes. 

From these considerations, we may construct the 


following empirical formula applicable in the range of 


ip Ste & | 
pen 7 (1030-F);1000> 8 >950 
H= oft+n coo(e-a)} \* te 
Wa 11.43 5; —&<950 


where H is the wave height in meter, RL the central pres- 


sure of typhoon in millibar, Y the distance from the 


storm center in nautical mile and 9 tne angle measured 


(C) C > 40 knfr 


clockwise from a line showing the direction of movement 
Fig-2 Mean distribution maps 
of storm center. € means the distance from the storm of wave height (in meter) due 


R to typhoons whose central pres- 
center to the region where the highest waves are raised, sures are 950~980 mb in either 


; ; case but whose speeds (C) differ. 
and a, n, a and % are constants which depend on the 
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speed of storm, though we have a=1.02, n=0.20, and A&=135°on the average. A calcula- 
tion diagram based on this forrula was already published. In the range of Y<€ , suffi- 
cient knowledge to construct an empirical formula has not been obtained yet, but qualita- 
tive discussion will be made in section 5. 

Moreover, it must be kept in mind that, out of typhoon too, we can generally find 


swells with greater size and length in the rear than in the front as well as within the 


' storm being contrary to what is usually believed. 


3. Distribution of wave period 


somewhat smaller than the observed one in shallow water 
near a land. Moreover, the wave period in the front semi- 


circle is likely to decrease as the speed of storm in- 


The distribution of period of prevailing waves is re- 
presented in Fig-3. Waves with great length are more fre- 
quently raised in the rear semi-circle than in the front 
semi-circle, and the wave period decreases on the whole as 
the distance from the storm center increases and as the 
central pressure of storm becomes high. The wave period 


frequently observed in typhoon is 7~11 sec, and this is 


creases. Fig-3 Mean distribution map 


However, it should be inpressed that actual seas consti- 


of wave period (in sec) due 
In conclusion, we may say that the distribution of to a typhoon whose central 
pressure is 940~960 mb. 
wave period is generally similar to that of wave height, 


although the former is more complicated and more confused. 


It we neglect the difference of wave period among quad- 


‘rants and the effect of speed of typhoon, the distribu- 


tion of period of prevailing waves is given apprcximately 


BR i fa 


ty a Ber OES ok 
T = 485-35 ~- 250 


tute complicated spectra including waves with various 


periods as shown in section 7. 


4. Distribution of wave steepness and wave age 


The wave is the steepest in the vicinity of the 


storm center with the steepness more than 1/20, and it 


becomes more gentle as the distance from the center in- Fig-4+ Distribution map of 
wave steepness in typhoon 
ereases (Fig-4+). Moreover, the wave is steeper in the whose central pressure is 


940~960 mb. 
rear than in the front. The steepness observed frequent- 


ly in typhoon is 1/15~1/30, and the steepness more than 
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1/10 is very rare. 

The wave age decreases usually as the distance from the storm center diminishes. For 
example, in typhoon with 960.980 mb the average of wave age is 0.60, 0.81 and 0.95 at the 
places whose distance from the center is 0~100, 100~*200 and 200~300 sea miles, respec- 
tively. Any wave whose age is less than 0.4 or longer than 1.4 (=6m) was not reported as 
a prevailing wave in question. The wave steepness for a given wave age is generally smal- 


ler than the value given semi-theoretically by Sverdrup and Munk (1947). 


5. State of the sea in the eye of a tropical cyclone 

It has been usually believed that in the eye of typhoon or hurricane pyramidal waves 
are raised because waves penetrate into this place from all directions and collide each 
other. On the contrary, Ootani (1941) suggested that pyramidal waves would not be raised 
in the eye, but they occur just behind the storm center. Then, the author collected data 
showing the state of the sea in the eye observed from ships or planes. From them it could 
be recognized that the sea in the eye is less confused when the typhoon moves slowly, and 
becomes more confused according as it moves faster. For instance, the sea in the eye was 
"moderate" when C (speed of storm) =7 km/hr, "confused" when C = 20 km/hr and "very con- 
fused" when C = 36 km/hr. 

The cause of this interesting phenomenon can be made clear by considering geometrical 
conditions necessary for the penetration of developed swells into the eye and also by con- 
sidering the effect of wind on these swells. Thus, the sea condition in the eye of tropi- 
cal cyclone is expected as follows: 

less confused when C<V, 

very confused when C~V, 

rather confused when C >V, 
where V is the group velocity of wave (usually 20~40 km/hr). Particularly, in the rear 
region not so distant from the storm center, "mountaineous and pyramidal waves" are apt 
to be raised when the storm progresses with a speed larger than V, because the remarkable 
swell originated in the front semi-circle collides at great angles with the wind wave gen-~ 


erated in the rear (confer Fig-6). 


6. Direction of wind waves and swell 

In relation to the presumption of the location of tropical cyclone, it has been dis- 
cussed a little which direction wind waves and swell in it take, but the discrepacies of 
opinions are found. Actually, the distribution of wave-direction is not so simple as 
stated formerly, but depends on the individualities of storm such as the speed, the track 
and the stage of development, etc. 


Fig-5 represents a typhoon which has but just developed rapidly. We recognize that 
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rig-5 Composite map of wind, 

wind waves, swell, period and 

height of prevailing waves in 

Tyohoon Jane of Sept. 2~3, 

1950. 
i waves (or wind waves of swell type) prevail 
a the same direction as that of wind, since 
lls have not grown sufficiently yet. Fig-6 shows 
ase when the typhoon moved straight forward with 
speed of 50 km/hr which is greater than that of 
il. In the front of center wind waves due to 
is at that place prevail, while in the rear of 
ter conspicuous swells which have been left be- 
i coexist with wind waves taking almost right 
les between them. When the storm advances slowly 
7-7), swells with direction different from that 
vind waves are seen in the front of center as 
' as in the rear. When the storm is changing its 
*se the distribution of wave direction becomes 


Jlicated, and when the storm is located near a 


\y 4 
28.03; 100,70 SS 


2809;;10.0.60 
28,1270, 5.0 
CS a ee 
Nautical Mer 
FigeS Composite map of wind and 


ocesn waves due to Typhoon 
Patricia of Oct. 25~28, 1949. 


24’ 


heel Vi 


Fig-7 Composite map of wind ana 
ocean waves due to typhoon of 
Oct. 1946. Small arrows with 
solid or broken lines indicate 
the direction of wind or swell, 
respectively. Numerals in the 
front or rear of the arrows 

give the scale of swell and the 
speed of wind (in knot), respec- 
tively. 


i the influence of being of land on the direction of waves is also recognized. 


Further, observed values of angle between wind and swell can be explained quantitatively 


‘espective cases from a theoretical consideration. 


When the storm moves slowly the swell 


th should appear in whole area deviates from the wind to the right (in the Northern Hemi- 


re) as pointed out first by Tannehill (1936). 
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As the speed of storm increases the 
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deviation of swell from wind becomes small in the right section ot storm and large in the 
left section. Particularly, when the speed of storm exceeds tnat of swell, we can find 

a swell with almost an opposite direction from the wind in the left rear quadrant, and in 
the right rear quadrant we can find a swell travelling to the left of wind contrary to 
Tannehill's consideration. 


In the case of typhoon associated with a frontal system confused seas and swell with 


unusual directions are likely to be raised in the adjacent seas of the front. For example, 


in Fig-7 we observe considerably strong swells deviating from the wind to the left on in- 
tense warm front which lies to the left front of center, and confused seas are found to the 


right rear of center where the weak cold front is located. 


7. Wave energy spectrum in typhoon 

Since it is difficult under present conditions to obtain wave spectra in tropical 
cyclone on the open sea, the author intends to investigate the distribution of them on the 
basis of wave energy spectra that Ijima (1957) constructed by analysing wave records at 
the Port of Onahama. Fig-8a indicates the relative position of the port. to the storm 
center at each observational time when the energy spectrum shown in Fig-8b and ¢ could be 
traced. In this case, we may say that the variation of wave spectrum with time at the port 
represents qualitatively the distribution of it in the right semi-circle of typhoon on the 
open sea, because the storm passed west ot Onahama and the speed of it is very great (60~ 
90 km/nr). 

In the right front quadrant of typhoon with great velocity, as it becomes near to the 
storm center, the peak of the spectrum moves rapidly towards the long period band result- 
ing in the increase of energy in this band, while the energy of short period band remains 


constant. it just agrees with results given by the tneory of wave generation. In the 


420 ESO Do ection of Movemext WS 0/0 Ws 020 fisey 
the typhoex. centre 
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Peoeee ae ec the time change of wave energy spectrum 
at the Port of Onahama during the passage of Typhoo 
Tess from Sept. 25 to 26, 1953. ; i es 
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.r region of the roughest waves, as the distance from the center becomes large, the 
eT gy decreases gradually at every period keeping the period of peak in wave spectrum 
be constant because swells are prevailing (see Fig-6). 

Moreover, on the open sea, it is also expected that the wave spectrum in the left 
ni-circle does not differ seriously from that in the right semi-circle. But, when a 
ohoon moves slowly the change of wave spectum shall be different from the above- 
itioned case. Consider, for example, a case when the storm is at rest and the distri- 
tion of wind is symmetrical about its center. The spectrum is to be naturally just the 


ne for all quadrants, where wind waves and swell coexist. 


Energy of waves raises in typhoon 

The total energy of waves amounts to about 100 age in a strong typhoon and less than 
? erg in a weak one. It corresponds tc about ten percent of the kinetic energy of the 
9hoon itself. 

If we assume that the state of the sea is steady and the dissipation of waves due to 
iction is neglected in accordance with Sverdrup and Munk (1947), the energy of waves 
nerated per unit time by wind in the storm must be equal to the flow of energy per unit 
ne due to swells which go out from the storm. From such consideration, it is shown 
at several percent of the wind energy dissipated by surface friction goes toward increas- 
> and maintaining the energy of waves, and this agrees nearly with what was obtained for 


snificant waves by Sverdrup and Munk from entirely different considerations. 


Mean relation between wind velocity and wave height 

When a certain velocity of wind (W, m/s) is observed in typhoon, the average of wave 
cghts in all cases (H, m), to be observed at the same time, is given by an empirical 
emula 


H=0.31W 


shown in Fig-9. Such linear equation is similar to those given by Cornish (1934) and 


ners, but our coefficient 0.31 is naturally smal- fe 2 eet y aie __ 42:Wlbe/s) 
> than those in the latter, since they are related 


the highest waves. In Fig-9, we recognize that 


> wave height in the Atlantic given by Roll (1954) 


almost always smaller a little than that in ty- 


yon though they are parallel, anc also that the 


sbable wave height (recommended by the World Me- 


yrological Organization) is smaller when the wind 


Fig-9 Mean relation between 


weak and larger when the wind is strong than the the wind velocity and the wave 
height in typhoons. 


ight in the typhoon and in the Atlantic. 
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Besides, when the wind velocity is increasing continuously the wave height is usually 
lower than that calculated by this formula, and on the contrary, when the wind velocity is 


decreasing continuously the wave height is usually higher than the calculated one. This 


can be interpreted by the fact that in the former case the weather ship is located in the f. 
semi-circle of storm and in the latter case in the rear semi-circle, therefcre, the waves 
are lower in the former than in the latter ("ig-1). The mean relation between wind velocit; 


and wave period is omitted here. 
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THE POSSIBLE HIGHEST GRAVITY WAVES 
IN SHALLOW WATER 


Tsutomu Kishi, Dr. Engineering 


Assistant Prof. of Hydraulics, Hokkaido University 


: When we design seashore structures, for examples, groins, breakwaters and sea 
c., it is often essential to consider the possible highest clapotis and progressive 
shallow water. 
present paper is an abstract of two papers published by the author in Japanese dur- 
and 1955, and consists of two parts. In the first part the clapotis of finite height 
he second the progressive wave of finite height are respectively discussed. 
criterion of the highest clapotis is based on the vertical acceleration of water 
» while that of the progressive wave is based on the horizontal velocity of water 


e 


possible highest clapotis in shallow water. 
ution of the clapotis of finite height in shallow water. 
x-axis be horizontal and y-axis vertically upwards. Velocity potential ¥ and wave 
(x, y,t)=0 are assumed to be expanded in the following infinite series : 
P =D anlt) cook nk(y +d ) cosnk x (1) 


n=0 


oe ( 
VOY t)= YZ Qo(tt ® An ltesonkx =o, zt 


bove expressions, a,(t) and @,(t) are unknown functions of time t and k is the wave 


Equation (1) must satisfy the condition of continuity. Namely, 


a7Y ay 
pee tneng ye TNCs, (3) 


the wall surface the horizontal velocity u is always zero so that 


ae : (4) 
ox (9) ate -1Oee LOM allt <. 


SS 
" 
( 


resent problem, the water is of uniform depth d. And the vertical velocity will 


at y = - d, i.e. 


ae eT AtaymO forvald, t. (5) 


og 
itions of eqs. (4) and (5) are satisfied by the potential assumed as eq. (Ge 


m the condition of pressure p=O at the free surface, one has 


-9- 


- OF _ / OP es g 2 = =O 
dais ae Ee + (3f) 0 pag ay (Cg Boe 
At the free surface a fluid element must move so that its velocity component normal to ty 


surface is the same as the normal velocity of the surface itself. Hence the function 


must satisfy the condition 


ne ee on hg, t) 0 (7) 


Therefore, the problem is reduced to the determination of the unkmown functions a>(¢t) a 
Q»(t) under the conditions given by eqs. (6) and (7). 
By substituting eqs. (1) and (2) into eqs. (6) and (7) and equating the sum of the 

respective terms of cagn/px to zero, following equations are obtained: 
from condition (6), 

£905 + be + h(x, a, airh kd + £4, 05 Q, Cooh Rd + 222 Ay 201 

24% aE Cat aot ge 71 bo Xi + 2h, A, anh 2RkA 

+ 2Rdz & A, cosh 2kd. + Fs Qy Ainh3 kd +£RX,Q,Q; Cook 3RA ) 

- 5G: Cabh 2kd + RAP QAo anh 2kd )=0 

- 9a, ta, Coohkd+ BA, Q sinkkd + #543 Coohkd+ Bx, Q, anh ka 

+ BH, QoQ, CORA + R&A, AUN 2KA + Rol Qos Coa 2RA-k’adpo)\ CodhIkA = O (8) 

- 92 +02 Lod 2kd + RM sinh 2RA tH HQ) 2irh Rd + Hol Ao ds COSA 


+ ja; auch Rd + R73 Ay Qs Bonk +H 02 C09 2KA t PREY Ao ee SRAt R= 0 


= G83 +A; Cod 3kA t BRO 32o atnh3kd + $ R72, cooh 3kaA 


+ £8, a, anh kat 4 Aods Coon kd thes Q, sinh 2kA 
+ K?, 2,8, Cah 2kd + Ror0l, Coon Rd =O 


and from condition (7), 
La, ~- & (x, 2, cok hdrby,a,2, 4ink RA+ 9K;3A, Cork 3hA 
+ 22 kK3 0,4, Anh FRA +4X2Qr Cork 2RA +4hd2QyQ, sink 2A.) =0 


by tha sink hd + Bia, ap coohhd ~hi(o/Qx Coch kd hot, Aol, aink hed 


+ &,Q, Cork 2hd +RhA,Qaamnk 2rd) =0 
QA2t+2hk, Ankh 2hd +RQ, cosh hd + 2h, Qo cook 2hd 


+H (Bick AoA, rink hd + dol, Q, ooh 2kd + Aol ls ty sink apd (9) 


_ lo- 
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+ 903, oar sikd+ 27 pry, 2, @, Ankh, Jka )=0 


Ast 3RAK; nk IRA + 2 holy Q cook Id +f*(0,A2 Cork hd 


a 
| +2 Qa, tink Rd + ola Q, Caak Lhd + fel, Qo @, anh Ded) =O 


) oS Sa a ee et mee 


cessive calculations give the following results: 
Keer C= Az oiain? + 2 Cosh 2hda_ 2 
29C-> ot +£AG? Gate COs o t 


A= = Cooo t 


AG 
Ramh kad 


Eee: 
i 2aink Akd (2tanh hl -tanh 2hA) 


Qo =C 
IO) 


= A sm o-t 
74 4 dtn* hd 
3kA? tank kad 


ee SORA lore RR 
 Qaink 2hd (2tankhd- lankAZRdad ) 
= Pk tank kd (/--4£ tanhhdt Camb hd + age 


Where A= const: C = conat., 


2) 


Cag 2et 


> physical interpretation of the first constant A is easily obtained by considering 


-itude of this wave. 


Vmax A +f tank ka { 


fe Z 
ainhihd( 2tanrhhd- tanh 2rd i. 


=A + BE tank kd | 
ie AimhZka Cs larh 2hkA) 
means the half of the wave height of clapotis. Second eee C indicates the 


alal) 


Umi 


of the height of the still water level above x-axis. Referring to the Sainflou's 
1, the author assumed the value of C as 

C=fkA? coth ka, Gen 
> profile given by eq, (10) has been examined by the author to show a good agreement 
> experimental results reported by F. Suquet and A. Walletel 
> highest clapotis. 
the perfect liquid cannot bear the tensile stress, the condition of op/ay < 


”e satisfied at the free surface. The equation of motion related to y-axis is 


y {Fo ‘ 
eet sy ta a ca are (13) 
=O at te free ae when the wave is critical, the above equation becomes 
av 2 aay He Zi 
ott Poy ae 


) have proposed the following relation as the criterion ofj 


W. G. Penny and A. T. Price” 
the highest clapotis in deep water: 
av = (14) 
Epa res 
And the maximum height calculated by eq. (1+) is 


RA = O29 2 

In the case of a shallow water, the maximum height will decrease as the water depth 
decreases. Then the generalized criterion of the highest clapotis in any water depth ca 
be expressed as 


20 5 Ac Gen oe & is a coefficient and S&S / (15) 
The coefficient x in eq. (15) must satisfy the following condition: . 


1. Lem kA =0 


Pe, Lome, RA = OS 92. 
The author assumes the coefficient Xx in the following form: 
ad = tank? ka. (16) 
The relation (16) is one of the simplest relation that satisfies the above two condition 
Let the origin be on the wave crest and then 9Y/g9t becomes from eqs. (10) and (12) a 
v i iy 
ae tee = — gGkAtanhed{) + hh coth kd + o(Rg2)t- 7) 
Substitution of the above result into eq. (15) leads to 
2 
(RA coth red)” + (kA coth hd )-/ =0 + 0 (RA?) (18) 
As the expression (17) is applicable only in the order of magnitude o(r7a*), eq. (18) 
should be replaced by 
2 
(kA coth kd + (KA cethhd)-~ = 0 (163 


Where is a coefficient nearly equal to unity. 


By determining the value of @ in such a way that the solution of eq. (19) satisfies the 


’ 


condition 2, i.e., Ax. kA = 0572 ,68 = 0.940 and the solution 
kActhAkd = 0592 (20) 
are obtained. The relation (20) is the general criterion of the highest clapotis in any 
depth. 
In the case of a practical application, it is convenient to express d in relation 


(20) by the water depth H. When the wave is critical the relation between d and H becon 


a=H-F= H-ZkRA cothhkd = H-A/4sA 


(21) 
Substitution of the relation (21) into the relation (20) gives the relation 
kA = (cow? kH t 0.350 Cesech*kH) - Coth RH ; 
0.296 cadech*? RH : k22) 


The relation (22) is compared with the experimental results of F. Suquet and A wallee 


hat Boke, “Dkr. 
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© Calculated value 
e Measured value just before wave breaking 
+ Measured value just after wave breaking 
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rding to Fig. 1, the criterion of the highest clapotis given by the relation(20) or 


seems to be applicable. 


The possible highest progressive wave in shallow water. 

Solution of a progressive wave of finite height in shallow water. 

In the author's method, the progressive wave with celerity c is superposed with a 
orm flow of equal velocity in the opposite direction by following the Rayleigh's method. 
resulting motion becomes a steady flow with a sinuous surface. 

Let x-axis be horizontal being positive direction in the superposed flow and y-axis be 
ical upwards. The velocity potential f and stream function ¥ of the flow under consid- 


ion can be expressed as follows: 


fp =-CXt+CHRcopk(Aty )ank< , (1) 
b= -CY+CB 2h k(dt+y) Cookx 
where k = wave number = 20//, 


L = wave length — 


c = wave celerity 


Stream lines Y =O and &% = cd indicate the free surface and the bottom respectively. 
profile of the free surface given by 5S. eee te 
yh =$hacothhkd + acoohx +Lharceth ha cos 2X 
+ $h?Q? (Coth kd + $) cov3 wx 


where 
displacement from the axis y=0 


uy 


Bamk kh {i+ 4 kB? (/ +f aemh? ht)} (3) 


larly, the mean level y of the strean line 2 =ecqd'! (a'< a) and profile of the stream line 


a 


ziven by eqs. (4) and (5), respectively. 
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y= -#'+$ yc hk (d-a’) () 


b'= +r? coth kh (d-d’) +y cookx+ SRy?eobk(d- a’) cod 2Rx 


+ fer | cob*e(d- a’) + $j covdkX torres (5) 

where : 
es anh k (a =-d’) 4 
\'* paca vepdenhohiias wa SS 


These relations give the momentum ccentained in the domain of unit length bounded by 


the stream line 4% = O and &% = cd'. The momentum 77s50f the steady flow is given by 


¥=0 
Ms Rr Ud = is og, = a 
$= rae rare bcd’ OY dy = fea Ce 
where f density of the fluid — 
While, the momentum of the superposed uniform flow of velocity ¢ contained in the same 


domain is given by 
fC 44 (3) 


Consequently, the momentum 77y/of the progressive wave with the propagation velocity c is 


the difference of the above two expressions, that is 
’ ma! = fC(A-Y) CodkRX + Pek { arcoth kd -Ycoth aC d- a’) cook X 
++ pcr?| @( cthka + $) 7? 4 coth?R( d- L' + Zh | caa?ext™ ak: 
ele ra 2 
By putting d=d' in eq. (9) and by integrating over one wave length, the momentum of the 
progressive wave per one wave per one length can be obtained: 
M= pcka?ahrd = 2S CQ? cohkd . (20) 
Horizontal velocity(U(g/along the stream line ¢=Cd' can be calculated by differentiating 
(M/P) with respecting to d'. 
Ud’ = 7 Sak = Cak Gash EG 42 cookx + Ca%p? oes A c0a kx 
+ecae me Wee a dh cooh*k(d - a ee z aes LO Ts 16612.) 


tinh? -#a 
when d! = O is substituted into eq. (11), the horizontal velocity along the bees surface 


Uo becomes 
Uo =Cak cothkd cookx + Ca7h( lh *k d+ /) coo? RE + 4 CAR ?( coth® ed + J coth kd) Cood Xt 
The value of (/p takes its maximum at the wave crest x=0, i.e., (12) 


Uomax = Chak Coth kd + Q°k?( coth* kd Tr) 


+ Zak? (coth kd + J coh kd)+ >. J a8 
(2) The highest progressive wave 
A wave attains its possible highest limit when the horizontal velocity of the fluid 
particle at the crest becomes equal to the wave celerity. The above fact gives the fol- 
lowing criterion of the highest progressive wave: 
Ak cothkd+ ap coh? kdt+/)+Farri(cth’rdtfecthkd)=/ (14) 


In order to examine the above relation by exveriments, physical meanings of the two 
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ameters "a" and "d" involved in eq. (14) must first be revealed. From eq. (2), it is 
ily seen that “a and “d" are respectively related to the half wave height A and water 
th H by the following realtions: 
ait (15) 
H =ht+tRarcthrad 

The combination of the relations (14) and (15) gives the relation between steepness 
L and shallowness H/L of the highest progressive wave, which is shown in Fig. 2 by a 
id curve. The experimental rsults reported by F. Suquet and A. Wallet *are also 
ered in the same figure. 

The calculated values agree 
paratively well with the experi- 
tal values. Another examination 
eq. (14) is made by considering 


limiting case of deep water, 


by 
«, kKd=co . When the condition : Lasse 
= co is introduced in eq. (1+), ak= 0. iil SS WNL 
60 or 2A/L= 0.146 = 1/6.8 is (ea ae 
ulted. The above value of 24/L 


Ply coincides with the value of 


6) b 
given by Michell. 0 0.05 010 O15 0.20 025 0,30 0.35 0.40 045 0,50 055 


Fig. 2 
elusion 


In the present paper the author discussed the possitle highest gravity wave in shallow 
er. The ciriterions of the highest waves, clapotis and progressive wave, proposed by 
author agree fairly well with the experimental results. However, he still has some 
stions from the mathematical point of view because the convergences of his solutions 


not verified. 
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FIELD INVESTIGATION OF WAVE FORCES 
ON BREAKWATER 
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ntroduction 

Many theoretical studies and laboratory experiments on wave forces against breakwaters 

been made. However, it is usually very difficult to apply these results directly to 
natural sea because of the complexity of the actual phenomena. 

The present investigation is intended to make a comparison between experimental results 
ined in laboratory and the data observed in the natural sea, and to determine a practi- 
method for designing economical and safe breakwaters. The practical observations of 

dimensions, wave pressures and breakwater vibrations have been made at Haboro Harbor 
2 November, 1957. 

This paper presents an outline of the experimental equipments and some of the experi- 
al results obtained -- the classification of wave pressures in shape, frequency of 

type pressure, some considerations on the maximum wave pressure and vibrations of the 
kwater. 
sscription of Haboro Harbor 
Haboro is a small fishing harbor located 
ne north-west coast of Hokkaido facing on 
Japan sea. This harbor was designed and 
© as a river-mouth-harbor in early time, 
now, the course of the river has already 
changed artificially as will be seen in 
|. At present, the construction of new 
cwaters which are shown by the broken line 
ig-1, is going on under a new plan for the 


»se of preventing-waves and drifting sands 


entering into the harbor. ee, 

This harbor has several convenient points | ioe 
investigating wave pressures: HABORG HARBOR pee ey! 
) It is easy to set measuring instruments Fig-l 
on the breakwater, because the new 
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breakwaters are mow being built. ] 
2) Since violent stormy waves come to this sea region during the period from November 
to February, many available data for investigations can be expected. 
3) Height of tidal wave is very small. | 
4) The direction of the breakwater which is used for the observation is perpendicular . 
to the direction of the main wave. ; 
5) As the foundation-mounds of the breakwaters are new, the up-lift pressure which 
acts against the bottom of the breakwater may be measured easily. 
6) The obtained results can be used directly for the practical works of breakwaters 
at this harbor. 
3. Equipments provided for the observations 
The principal instruments provided for the measurements of wave heights, wave pressures 
and vibrations of the breakwater are as follows: 
(1) Pressure type wave recorder 
This wave recorder was designed by Transportation Technical Research Institute, 
Japan. The principle and mechanism are similar to California Mark III wave recorder. This 
wave recorder was provided for measuring the offshore wave heights. 
(2) Transit type wave recorder 
This wave recorder was designed by G. Udai, one of the authors. The principle of 
this wave recorder belongs to the optical method. First, a float which has a clear mark 
is set at the expected position of the sea. The motion of the float, which moves closely 
with the waves, is traced by the observer with the transit which is set in the observation 
room built at the higher place. The motion of the transit is transmitted to a recording 
chart through a magnifying lever mechanism, and the wave motion is recorded with pen on 
a chart. This instrument has some excellent points. The principle of this instrument is 
very simple, and it seldom be damaged. And, its cost is very low. But it can not be 
used when the weather is so bad, because the mark of the float can not be clearly seen. 
This wave recorder was also provided to measure 
the offshore wave heights. This apparatus is 
shown in Photo-l. 
(3) Step type wave recorder 
This wave recorder was used to measure 
the heights of the waves at the front position 
of the breakwater. This recorder was designed 
by Y. Muraki, one of the authors. Similar 
type wave recorder which has been used up to 


the present involves an inconvenient point -- 
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shough relatively clear steplike wave records can be obtained on a chart when the sea 
-e rises, the records show smooth curved without showing steplike shapes on account 

> electrical leakage through the sea water film of the electrodes (spark-plug) when 
2a surface falls. It means that the records accompany the time lag and do not show 
-tual wave motion. 

fr. Gerhrdt also noticed the same point in his report. He used the electrodes which 
mated with silicon grease in order to avoid the influence of this leakage, but it was 
‘ficient. 

Ine of the authors solved this point by using magnet-relays, which have .ow electrical 
pances, on the old type apparatus. The detailed explanation of this gage will be re- 
1 in another paper. The diagram of this circuit is shown in Fig-2 and an example of 
scords, which shows the very clear steplike shape from start to end, is shown in 

“26 


pa hl IN le Ip sea Ee RN IERIE be DET ES ee 

Jave pressure gage | 

SS ee | ELECTRODES| 
har 


ilectric Resistance Strain Gage SHO O er 


Jave pressure gages are used in 
investigation. These were de- 


i by Transportation Technical 


‘ch Institute, and made by Toyo 


| 
1e amplifying apparatus are | x Bi aS 
6V BATTERY 
in Photoes-3 and 4+, respective- | i 


ibrograph 


K.K. The pick-up of pressure 


mn order to investigate the mo- 
f£ the breakwater caused by wave 
res, and accelerograph is set 

. breakwater. At first, the 

s had an intention of measuring 
bration in displacement, but a 
le displacement vibrograph was 

und . Therefore, a preliminary 

ation of only the period of the 


ion is performed with an accel- 


ph in the present investigation. 


celerograph used is Electric Photo-2 


ance Strain Gage type which was 
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Photo-3 Photo-4+ 


STEP- WAVE-GAGE JOINT BOX & ACCELERQGRAPH 
PRESSURE ~ LEADING WIRE 
SUBMARINE-CABLE 


PRESSURE PICK-UP 


Fig-3 
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i. by Kyowa Musen K. K. Its proper frequency is 22 c/s and the capacity is + 1G. 

These measuring instruments are set on one breakwater as shown in Fig-3 and Photo-5. 
ur pick-ups of pressure are set on the front vertical wall and three pick-ups of pressure 
-e set on the bottom surface. The step wave gage, 6 m in length, is also set on the same 
‘ont wall and the accelerograph is fixed to the breakwater in the manhole. 


: There is a joint box in the manhole. Leading wires from each pick-up of pressure, 


sep wave gage and accelerograph are connected together in the box with a submarine-cable 
1ich comes from the observation room 300 m apart from the breakwater. 
, Types of wave pressure 
.) Classification of wave pressure in shape 

The practical field observations were carried out for about three months from November, 
957 to January, 1958. During this period, many stormy waves which involved wind waves and 
wells came to this coast. Therefore, many different kinds of wave pressure records, in 
agnitude and shape, were obtained. An example of their records is shown in Photo-6. Ac- 
ording to the wave data obtained in past several years, it is evident that the season in 
1ich our observations were made takes the largest frequency in the occurrence of violent 
aves through a year. Therefore, the results obtained in this observation seem to repre- 
2nt the characteristics of the wave pressures, and the largest value of the observed 
ressures also seems to show an approximate 
alue of the maximum wave pressure to be ex-= 
=cted at this harbor. 

According to the records obtained, the 


ave pressures (at S. W. L.) take various 


napes. Of course, this variation depends on 
1e characteristics of the wave, topography of 


1e sea bottom, tyne of the structure and oth- 


"Ss. However, when the variation of the 


1apes is considered, the ohtained records 


sem to include the shapes of almost all the 


nd of pressures which would occur at any 


sher sea obstructed by a vertical breakwater. 


The shapes of the obtained wave pressures 


lich act on the vertical wall at S. W. L. 


re classified as shown in Fig-4. Although 


is manner of classification is not always 


CLASSIFICATION OF WAVE PRESSURE SHAPE 


e generalized one, it is believed that all 


Fig-4 
wwe pressures in any other sea can also be 6 
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classified by this manner. 

As will be seen in Fig-4, the curves of any group of A, B and C are always smooth and 
simple, perhaps because they are caused by nonbreaking waves. Om the other hand, the curve 
of group D has a sharp angle, which is perhaps caused by breaking waves. There are some 
differences among the groups of A, B and C. In type A, the curves incline ahead, in type B 
they are symmetrical, and in type C they incline behind. The shape of type A is usually 
caused by wind waves, type C by swells and type B by both. The wave pressure of type D, 
which is caused by breaking waves, has been studied in laboratory experiments. However, 
it has rarely been studied in the natural sea. 

(2) Frequency of each type of wave pressures 

In order to discuss the stability of the breakwater, it is very important to know the 
intensity of the wave pressure, but, it is also very important to know the frequency of 
each type of wave pressures in the natural sea. 

The frequency of each type of wave pressures at Haboro Harbor is shown in Fig-5. The 


values of frequency are respectively 45 @ in type A, 


27 % in type B, 2+ % in type C and & 4 in type D. It 50 
is noticed that type D, which is caused by breaking ‘ 49) 
waves and considered to have particular characteristics, * a0 
shows very small value. The detail of the frequency is 5 20 
shown in Fig-6. Fig-7 shows some relations between wave ae 


characteristics and the frequency of type D. It means 
that wave pressures of type D occur more easily when the TYPE 
wave period is about 8 sec, and the intensities of the Fig-5 


wave pressures at that time are relatively large. From 
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e facts, it may be said that wave pressures of type D would 

ly appear when some conditions concerning the wave character- 

cs, structure and topography are superposed at the same time. 
The most typical pressure of shocking type observed is shown 


ig-8. The shape of this shocking pressure is almost same as 


WAVE PRESSURE ATS.W.L ton /nt 


one which is obtained in laboratory experiments. The time 
hock is 0.07 sec, the magnitude of the pressure is 11 ton/m?2 
ifle and 3.2 ton/m*in bourrage. The numbers of such shocking 
sures observed are only 2 or 3 in our data. Although the 
uency of this shocking pressure is very small, the pressure 


his type can not be neglected, because such pressure shows 


ACCELERATION x10°G 


lly the large value in magnitude. 


The fact that laboratory tests and field observations showed 
me shocking type of wave pressure indicates the validity of | 
bis 


similarity between the laboratory experiments and the natural 


03,00 A.m, 15, Dec, 1957 
omena. 


The maximum wave pressure at Haboro Harbor 

A great number of observed values of wave pressures are plot- aaa 
against the wave period in Fig-9. This shows a relation among the wave period, wave 
sure and wave type. According to this figure, it can be understood that the largest 


pressure which is obtained in this observation 


was caused by the breaking wave which takes 


the period of about 7 sec. From this result, 
1 

2 alas ge it may be predicted that the maximum wave 
F X SWELL 


pressure, which can occur at Haboro Harbor, 
will also be caused by such waves. 

In Fig-10, wave pressures are plotted a- 
gainst the rising time, which means the time 
from the start to the first peak of each wave 
pressure curve. The upper limit of these plot- 
ted values is represented by one curve, which 
shows that the wave pressure increases steeply 
as the rising time becomes short. From this, 


it is considered that the maximum wave pressure 


will be caused by the breaking wave type which 


WAVE PERIOD Sec 


takes very short time of shock, and the value 


Fig-9 of the maximum wave pressure may be determined 
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Fig-10 
‘ 


approximately by this empirical curve. From this figure, the value of the maximum pressure 


which can occur at this harbor can be read as 14~15 ton/m* 

,When this maximum value of wave pressure is taken, the wave height can be calculated — 
as about 10 m by Hiroi's formula. When the past data of waves and the actual topography 
are considered, the wave of 10 m in height can not occur at this harbor. On the other har 
according to Minikin's formula”, the wave height becomes to be about 1.3 m, which is too 
small for the actual condition of the sea at that time. From the consideration described 
above, it can be seen that both formulas are not suitable for the prototype breakwater at 
Haboro Harbor. 

6. Vibration of breakwaters 

In order to study the stability of the breakwater, it is necessary to know not only 
the wave pressure but also the vibration of the breakwater caused by the wave pressure. 
But inspite of its importance, the measurement of the breakwater vibration has rarely been 
made up to the present. A similar measurement was once made at Algiers Harbor about 25 
years ago. In that case, the three-components accelerograph of crystal type was used. Ac- 
cording to that results, a proper frequency of the breakwater was not cbserved within 1~ 
1000 c/sec, but the shake with the same period as the 
wave period was recognized. 

In our observation, an Electric Resistance Strain 


Gage type of Accelerograph was used. Its proper fre- 


FREQUENCY Number 


quency is 22 c/sec. The period of the breakwater vi- 


bration obtained in our observation is about 0.2 sec 


O15 0,20 = Q, 1.30 
as shown in Fig-ll. This result is different from the PERIOD Of VIBRATION Sec 


result obtained at Algiers Harbor mentioned above. Fig-11 
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By comparing the observed period with the proper frequency of the accelerograph, it 
an be seen that the obtained records represent the acceleration and not the displacement. 


| 
ne displacement is more essential for studying the motion of the breakwater. ~But it was 


ot obtained. Only the existance of the vibration with the period of about 0.2 sec was 
und in this observation. 

| The largest acceleration was caused by the largest wave pressure. The record is shown 
n Fig-8 with the largest pressure. 

From the observation, two qualitative relations between the wave pressure and the 
reakwater vibration were obtained. The one is concerning the phase and the other is con- 
erning the magnitude. 

The vibration occurs at the regular positions against the one wave pressure curve as 
hown in Figs-12 and 13. In typical case, it occurs at each position corresponding to the 
tart, midway between the start and the peak, every peaks and the end of the wave pressure 
urve. There is one interesting phenomenon. In the case of a shock=-type pressure, small 
ibrations of pressure almost always appear at the position of bourrage, corresponding to 
he vibrations of the breakwater, as shown in Fig-8. 

The magnitude of the acceleration is plotted against the wave pressure at S. W. L. in 
igs-14~17. From these figures, it can be seen that the upper limit of the acceleration 
akes different value in each type. In A and D type, the limiting value increases linearly 


ith increase in wave pressure. Especially in type D, this trend of increase is remarkable. 
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On the other hand, in B and C types, it increases 
little with increase in wave pressure. 
7. Summary 


In this paper, the authors discribed the pur- 


pose of this investigation, some explanations of - naa 


10 
WAVE PRESSURE AT SLL ton/mn* 


measuring equipments used, and some results ob- 
tained in this observation. The pricipal quali- Fig-17 

tative results are as follows: 

(1) Quite a wide variety of the shapes of wave pressure can occur, but they are classi- 

fied into four types which are named by the symbols of A, B, C and D, in this paper. Ac- 
cording to another way, they are classified into two types. The one is the nonbreaking 

wave type (steady wave type), which includes A, B, and C types, and the other is the break- 
ing wave type, which is named by type D in the present paper. 

(2) The values of frequency of each type A, B, C and D are 45, 27, 24, and 4 %, respectively. 
Although the breaking wave type D takes only 4 %, it should be noticed that the pressure of 
this type will easily occurs if some conditions are satisfied at the same time. 


(3) The typical shocking type of wave pressure which was obtained in laboratory tests was 
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30 Observed in the natural sea. 
) The maximum wave pressure at Haboro Harbor will be caused by the shocking type wave 
ich has the period of 7~9 sec, and its magnitude will be 14~15 ton/m? 
) The maximum pressure is of the shocking type pressure. However, its acting time is 
-y short. Therefore, in order to investigate the stability of the breakwater, further 
vestigations should be made under the consideration of its acting time «nd influences 
on the breakwater. 
) The observed period of the breakwater vibration which is caused by the wave pressure 
about 0.2 sec. 
) The vibration of the breakwater occurs at regular positions of one wave pressure 
PVE. 
) D-type wave pressure causes larger acceleration of the breakwater vibration than the 
ner types, and the value of its upper limit increases linearly with increase in wave 
essure. 

At present, the further observations are to be continued. The available results con- 
rning the distribution of the wave pressure and the relations among the wave pressures, 


ve dimensions and breakwater vibration will also be obtained. 
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Introduction 

When a sandy beach is subjected to attack for long hours by the uniform wave action 
efinite properties, the definite beach profile corresponding to the wave steepness will 
stablished. Since J. W. Johnson, bes expressed this equilibrium profile as a dimension- 
form, a lot of researches on this profile have been reported. This profile is not 
influenced by initial wave steepness, but also by the size of bed materials and its 
inal profile. Among these factors, Y. Iwagaki and ae pavarae ines recently reported 

t the effect of the size of bed materials. But the effect of the original profile is 
yet much clear, though it seems to be very important to make clear this effect when 
shange of beach profile is studied from the standpoint of the equilibrium slope method. 
nis paper the authors are to state some experimental results on the effects of the size 
ed materials and the original profile of beach upon the final equilibrium profile of 


mn and sand movement. 


Zquipment and procedure 

The experiment was carried out in a water tank of oreganic glass, stiffened with steel 
= works, 50 cm wide, 50 cm tall and 20 cm long. The methods of measurement of wave 
erties and collecting the sand transport are nearly the same as in the preceding Aupenee 


saxperiment consists of following two parts. 


rpiment 1. (Case 1) 
Tests were conducted on beaches carefully molded to two different initial slopes, 
ly 1:10 and 1:20. Each beach was molded from two kinds of bed materials, whose median 
sters were 0.2 mm and 0.8 mm, and whose specific gravities were 2.56 and 2.65 respec- 
ly. The wave properties submitted to the experiment were 80 = 210 cm in wave length 
sep water, 2.1 - 7.1 cm in wave height and 0.0104 - 0.0790 in initial wave steepenss 
Case 1-1. D=0.20 m, initial slope I,=1/10, {'=0.0107 = 0.0672 
Case 1-2. D=0.20 mm, initial slope Ip=1/20, 3;=0.0108 - 0.0719 
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Case 1-3. D=0.88 mm, initial slope I,1/10, §,=0-0109 - 0.0613 
Case 1-4. D=0.88 mm, initial slope Ip=1/20, §R=0.0104 - 0.0789 
The shape of beach profile was measured at every one hour in the case of 1:10 slope, 
and at every two hours in the case of 1:20 slope. The duration required for obtaining 
the equilibrium profile was ranged from 4 hours to 10 hours in the case of 1:10 slope and 


from 8 hours to 14 hours in the case of 1:20 slope. 


Experiment 2. (Case 2) 

In this experiment, the beach of 1:10 slope was molded, in the first place, from the 
sand whose median diameter was 0.2 mm, and then the beach was attacked by uniform wave 
action of a wave steepness on After the equilibrium profile was established by this wave 
action, this beach was attacked by uniform wave action of another steepness §;. New beach 
profile was measured at every one hour until the equilibrium profile was obtained. When 
the equilibrium profile was established, the shape of beach profile and the amount of 


change of shoreline were measured. The range of this experiment is as follows: 


» “Case 2-1.9° Initial slope’ 1:10, ""},=0.067, » H,=5.33,, Sy=0-0104 - 0.0456 
Case 2-2. ” $,=0.011, H,=3.55, }2=0.0157 - 0.0666 
Case 2-3. j §,=0-017, Hy =3.98, 5.70-0101 - 0.0647 
Case 2-4. s $.=0-022, Hi=3-97, },=0.0086 - 0.0598 
Case 2-5. & 5,70+032, H\=5.02, $=0.0090 - 0.0698 


3. The effect of the size of bed materials and the initial slope 

upon the equilibrium profile 

Fig. la, 1b, le, 1d show dimensionless expression of the equilibrium profiles for 
the different wave steepness to make clearthe effects of the size of bed materials and 
the initial slope of beach. The distribution of sand transport along the bed is also 
shown in this figure. From the figure we can see the difference in the location and 


the scale of longshore bar as the effect of the difference of the size of bed materials 
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cen 0.88 mm sand and 0.2 mm sand. Namely in the case of 0.2 mm sand, as the wave 

pness becomes steeper the suspended sand deposits in the offshore far from the break- 
point and forms large longshore bar. Local erosion and small bar are also seen near 
foam line. On the other hand, in the case of 0.88 mm sand, the range of sand movement 
imited in the parts between the shoreline and the breaking point. In the case of low 
steepness the amount of sand transport is concentrated near the shoreline. In the 

of large wave steepness, the distribution of sand transport of 0.2 mm sand is flat 

all range of movement, but that of 0.88 mm sand has two peaks, one is on the shore- 
and the other on the bar. These results agree with those of Iwagaki's experiment. 

In the next place let us consider the effect of initial slope. In the case of 0.88 mm 
. the effect cannot be clearly conceived as the range of sand movement is limited be-~ 

1 the breaking point and the shoreline. In the case of 0.2 mm sand the beach profiles 
sen the shoreline and the breaking point are nearly similar, but the profiles in the 
10re show some difference according to initial slope. 
Fig. 2 shows the relation between the initial wave 
sness and the dimensionless expression of the amount 
and transport at the shoreline. From the figure we 
see the sand transport decreases rapidly as do becomes 


iffect of the original beach profiles upon 


she equilibrium profiles 


P 5 ein nama 

In this experiment the beach of 1:10 slope was molded o% i5 
2 As 4+ —-+— 
the sand whose median diameter was 0.2 mm and then | 
. | o jeod 

Yeach was attacked by uniform wave action of a wave Diss ss : 7 > 
mess >,° After the equilibrium beach profile was 
lished by this wave action, this beach profile was ules 2 
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Fig. 3a fig. 3b 


established by this wave action, this beach profile 
was attacked by uniform wave action of another 
steepness Sah Thus, the new equilibrium beach pro- 
file was established. Fig. 3a, 3b, 3c show these 
profiles. 

In the Fig. 3a the original steepness of u- 
niform wave action is 0.607 and the subsequent 
steepnesses are smaller than the original steepness, 


namely 0.01, 0,017, 0.023, 0.030 and 0.046. These 


are the case that the storm beach established by 


Fig e 3¢ 


large steepness b,is subjected to attack by wave 
action of smaller steepness than D1 - From the figure, it will be seen that the sand 

drifts deposit near the shoreline,the shoreline moves towards the offshore and the long- 
shore bar is eroded. The rate of this action is remarkable as the wave steepness becomes 
smaller. But, in the case of the wave steepness above 0.03, any changes in the shoreline 
or the foreshore are not seen. On the contrary, the surfzone is a little eroded and some 
deposits are found in the offshore. It seems that this comes from the movement of small 
bar made near the foam line of the original beach profile. We cannot see the remarkable 


change of the beach profile in this case. 


In the Fig. 3b, the original steepness of uniform wave action is 517 0-022 and this 
is such case that the beach lies between summer type and winter type. The subseqtient 
wave steepnesses are 0.010 and 0.013 which are smaller than the original steepness at and 
0.031 and 0.055 which are larger than the original steepness Dis In the case bes >» we 
can see some deposits in the foreshore, movement of the shoreline towards the offshore 
and the erosion at the foam line and offshore part. When J, is above 0.03 the erosion near 


the shoreline and the formation of bars in the offshore become remarkable. 
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Fig. 3c is the case that the beach established by the original wave steepness 5 =0.01 
subjected to attack by waves of larger steepness D2 than Se In this case we can see 
e erosion near the shoreline and the formation of bars in the foamline and offshore. 

From these results, the relationships between the rate of change of shoreline M/Ho 
i the subsequent wave steepness Sz2are shown in the Fig. 4 for the various original wave 
eepness Sie From the figure we can easily understand the shoreline moves towards off- 
ore or foreshore whether $2551. If 5,is given, M/Hovaries directly proportional to J2 

the range of §2$ 0.03, but it becomes approximately independent of 5, when $270.04, 

Fig. 5 shows the effect of the original profile upon the equilibrium profile. 
aches with initially different original profiles are subjected to attack by uniform wave 
tion of nearly equal steepness. From this figure we can see in the case }17 $2 the orig- 
al profiles still remain in the offshore part from the breaking point. 

From these facts, it may be said that the equilibrium profile consists of two parts, 
nely, the part mainly established by wave steepness and the part subjected to the in- 
uence of both wave steepness and the original profile. Thus, it comes that we cannot 
timate the equilibrium profile from the equilibrium slope method alone without any know- 
ige about the original profile. 

Fig. 6 shows the relationship between wave run-up and wave steepness. In the figure 
ye run-up is expressed by dimensionless parameter R/H,and experimental results by R. P. 
yagé are also shown for the sake of comparison. From the figure we can see that wave 
1-up has the definite relation to the wave steepness, but it is independent of the orig- 
al profile. 

Fig. 7 shows the relationship between the slope of foreshore and the wave steepness. 


> figure includes also the results by Iwagaki. From the figure it may be clear that the 
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slope of foreshore decreases as jo becomes large and it is independent of the original 


prifile. 


5. Conclusion 

In this paper, the experimental results on the effects of grain size of bed materials 
and the original profile upon the equilibrium beach profile are reported and discussed. 
The following are the main conclusions. 

(a) The equilibrium profile consists of two parts, namely the part established by 
wave steepness alone and the part subjected to the influence of the original profile. 
Generally speaking, the onshore profile from the breaking point is mainly subjected by 
the wave steepness and the offshore profile from the breaking point is subjected by both 
the wave steepness and the original profile. 

(b) The wave run-up and the slope of foreshore change with the wave steepness and 


are independent of the original profile and grain size of bed materials. 
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Synopsis: The purpose of this study is to analyze theoretically and experimentally 
the motion of water particle near sea bed, especially in surf zone, in order to obtain 
the fundamental knowledge on the dynamics of sand drift. The current meter of unbonded 
gauge type is used for the experimental study. "Yurthermore, this study deals with the 
theory of net mass transport velocity of water particles in surf zone. Introducing the 
idea of bed load functions in alluvial channels, it is possible to explain the exnerimental 
results on equilibrium beach profile fairly well from the change of effective water parti- 


cle velocity. 


ie Introduction 

It is one of the basic essentials to clarify the dynamic process of beach erosion by 
wave action for protection of beach and for pertient design of shore structures. Little 
literatures, however, have been published on the surf and bed material movement caused by 
wave action because of the complication in mechanism, and therefore no remarkable study 
nas been done on this subject. 

“It is evident that the relationship between the characteristics of sand movement and 

the ones of waves, and the difference of the amount of sand movement depending on the 
above characteristics are of basic significance to reveal the beach process, and therefore 
the hydraulic behaviours described above has urged strongly the author's experimental 
project. The dynamical analysis for equilibrium bed profile and transport mechanisms of 
beach sand, however, have scarcely been analyzed for their complications. Since the changes 
of beach profile result from the movement of sand drift, it can be declared that this 
problem is solved only when the mechanics of movement and dynamical analysis on equilibrium 
bed profile are established. 

There have not been remarkable studies on the mechanics of sand drift, but recently, 


Ippen and paaicen published several reports on the sorting of bed material, and the 
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study on movement of water particles produced by the motion of oscillating plane is in 

its progress by Kalkanis”, In these studies, it is revealed that the gravity force and 

the dynamic forces like drag, lift and so on produced by the progression of waves take 
main parts as influential factors on the sediment particle. Furthermore, making the as- 
sumption taht the velocity change of water particle by wave action is sinusoidal at all 
points, they conclude that this velocity change determines the movement of sediment par- 
ticles. Extrapolation of these studies for the flow on sloping beaches has not been 

made. It is not possible to analyze the mechanics of sand drift dynamically unless actual 
behaviours of fluid motion by wave are made clear, and consequently the author measured the 
movements of water particles at following five points on a model shore with smooth sloping 


bed by the use of the current meter of unbounded type” 


1) Offshore on horizontal bottom, 2) Inshore, 3) Near breaking point, 4) Near 

plunging point, 5) Near shore line. 

The measurement of spectra in velocity profiles by the oscillatory motion of wave 
were conducted at these points, and especially the following articles were investigated 
in details 

(1) Changes in maximum particle velocity in the direction of depth: The data ob- 
tained will be useful to get the characteristics of wave and furthermore its character- 
istics, which used to be studied only in terms of its shape and celerity, will become 
more definite. 

(2) Variations in maximumpparticle velocity along sea bottom: In this test, the 
transformation of incoming waves in their hydraulic characteristics is made, and the in- 
fluence of wave transformation on the sorting action of sand is estimated. 

(3) Change of water particle velocity in the course of time: With the sue of ex- 
perimental data obtained, the time variation of particle velocity in classical theory of 
waves, which is still popular for engineering use, is investigated, and furthermore the 
contribution to analysis of sand drifts will be possible. 

Using the spectra of velocity variation in wave flows the mass transport velocity of 
water, which will be largely related to the transport velocity of sediment particles veri- 
fied by Ippen and Eagleson, is obtained, and the back currents near breaking point revealed 
from variations in mass transport velocity along sea bottom. With changes in effective water 
particle velocity along sea bottom, the authors will explain the fact, approved by several 
researchers, that the sediment transport depends on the steepness of deep water waves and 
the grain size of bed material, assuming that the amount of sand drift can be expressed 
by the same expression as in the bed load in river. Furthermore, the authors intend to 


take into some considerations of beach process based on the idea of critical steepness of 
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6) 
‘es which has been proposed by S. Hayami. 


Apparatus for experiments and procedures 
Experimental flume used for tests. 

The pictorial details of the 

SS walled flume at the Ujigawa 


raulic Laboratory, Disaster Pre- 


ition Research Institute, Kyoto 


versity, are shown in Fig. 1. \ Pheumatic 
ave Generator 


consists of steel framed channel 


By Pass 


21.6 m long, 0.75 m wide and Valve 


2m high, and the pheumatic 


‘Ssure pump to generate the wave 


equiped. Steel bed of 5 m long Pneumatic Wave Generale 


0.73 m wide with variable slope 
placed at one end of the flume, Fig. 1 Details of wave tank and wave generator 
is forming a model beach. In this experiment, the slope of this steel bed was set at 
06', A by-pass valve is used for the pneumatic wave generator so that the pressure is 
othly adjusted, and thus waves generated in the flume varies their period from 0.7 to 
/ sec. The speed of revolution of the valve is regulated continuously by the Byer gov- 
rv of 1/21 motor and the wave height can also be adjusted by the by-pass valve. This 
erator is capable of producing the maximum wave height of 0.2 m. 

Current meter. 

The measurement of velocity variation in wave motion was made by the current meter 
eh transformed the difference of dynamic to static pressures into the electrical dif- 
ence of resistance. The actual velocity in wave motion is calculated by the output 
tage indicated in the meter. The capacitance gauge current meter has already been made 
M.I.T., whereas it seems to be unstable for measurements. Therefore, the meter is 
igned to transform the pressure difference into the change of electric resistance through 
_ strain gauge, The velocity changes are thus recorded in forms of continuous changes of 
ctrical resistance by the use of amplifiers and pen-*oscillograph. The principles and 
uctures of the meter is explained in next paragraph, together with its characteristics. 

This instrument transmits the strain of thin phosphor-bronze plate to unbonded type 
ge through a rod set at the thin plate and the resulting strain is transformed into 
_change in electrical resistance of unbonded gauge. The strain at the center of thin 
te which is subjected to pressure is expressed by the following equation which is an ex- 
ssion for deflection at the center of flat circle plate fixed to the boundary under 


uniformly distributed load (AP pressure difference in this case). 
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Sma =3-Y) 4P-V/16Et”, oi Mehaomts. etelitema Gls) 


where t: thickness of thin plate (mm), r: its radius, E: Youngs' modulus,)/: Poissons 
ratio, (t= 0.3 mm, r= 28 m in this instrument). 

Fig. 2 shows the details of the current meter of unbonded gauge. The portion of 
instrument described by (6) is interchangeable to either slender pitot tubes (3 mm 


diameter) or streamline headed flat type pitot tubes with slots (1 mm). depending on 


Pitot tube 
connector 


Connecting 
rod 


Fig. 2 Details of current meter 
of unbonded gauge type 


its purposes. Holes (9) and are 


air outlet of the meter. 


200% 


In all experimental runs, the pitot 


tubes of streamline headed type indicated 


: 


good results for the measurements of AS : aa 
ey max strain 3400x10 


water particle velocity in offshore zone 


where the pressure difference was small, 50 mex load /609r 


and the slender pitot tubes demonstrated Natural freg. 480°. : 
x10 


better results in the measurement of 0 2000 Strain 4000 


6 


current velocity than the former tubes 

Fig. 3 Characteristic curve of U gauge 
in onshore zone from breaking point where of T-3 type 
the pressure difference was so large that the irregular vibration caused by the surf gave 
much disturbances to the instrument. The U gauge of T-3 type was selected for this in- 
strument because of its better linearity, sensitivity and stability than other types. 
Fig. 3 shows the characteristic curve of this gauge. The maximum input voltage is 5 V; 


in- and output resistances 128.1, and natural frequency 480 c/s. The pressure load- 


strain relationship is shown in Fig. 4, and the maximum upper limit for measurement is 
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and pressure 
Fig, 5 Calibration curve of 


current meter 


gr/em*. The strain at the check point of 10 gr/em’* (pressure difference) is 1490 x Tee 


/mm. 

The calibration in velocities obtained by this instrument was made by comparison with 
ose obtained by the normal pitot tubes in the uniform open channel flows before conduc= 
ng the experiment. 

Comparison of both velocities is indicated in Fig. 5, where the solid line indicates 
e perfect correlationship between two values. It is found that this instrument gives a 
ttle larger value than the normal pitot tubes and the difference is less than 2 %. 

) Recording system and accuracy of measurements. 

All records were taken by the amplifier system at Disaster Prevention Research Institute, 
oto University, and the system consists of pre-amplifier, D.C. amplifier and ink-oscillo- 
aph. The pre-amplifier contains an oscillator circuit, an amplifier and a phase discrimi- 
tor circuit, and amplifies 2,000 times the original input, while the D.C. amplifier ampli- 
es 400 times, so that the system gives 800,000 times amplification totally which means 118 

gains. The accuracy of measurements depends on sensitivity of the recording system and 
aracteristics of the U gauge, and the later has a linear relationship as in Fig. 3 even 

jer low load. Thus, the sensitivity of the recording system is comparatively reliable, and 
makes possible to measure flow velocity of 5 cm/sec within errors of less than 1 %, using 
2 above amplifier system that operates an ink-oscillograph with full scale of 20 m at 200 
input, if the possible reading on the recording paper is as small as 0.5 mm. The time lag 
tween the current meter and the recording system was scarcely noticed after many checks. 

> results obtained with this whole system were good. 

) Wave height meter, 


The water level meter of resistance type was used in this experiment. 
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(5) Wave used for experimental runs. 
Three kinds of experimental waves were chosen among a series of waves. The wave length 


were 4.46, 2.84 and 2.08 m and the heights 4.62, 6.02 and 6.81 cm when transformed to the 


deep water waves. 
(6) Observation points. 

Current velocities for above three waves were measured at points from @ to @) shown 
in Fig. 1. The point @ locates in the region where the shallow water wave (h/L= 1/10 - 
1/2. where L is the wave length at the water depth h) develops and Q), in the transitional 
region where h/L = 1/10, and QB) - @), near the breaking points, (5) - ‘@) » in the region 


where the up-and backrushes prevail after breaking. 


3. Analysis of experimental results and its consideration 
(1) Variations in maximum velocities of water particle in vertical direction. 

In classical treatment for engineering uses, wave heights, lengths, and celerities have 
been only subjected for theoretical and experimental studies on characteristics of waves in 
shallow waters and surf zone that might be of most important factors on the movement of sand 
drift. Although many studies on the variations in maximum velocities of water particles wit 
the change of wave characteristics in surf zone, based on the wave theory in shallow water, 
have been made, it is impossible to discuss the movement of sand drift unless the water 
particle movement by wave action is studied, since the movement of sand drift is essentially 
a dynamic process in behaviour. From this stand point of view, the authors have made an ex- 
perimental research on the changes in maximum horizontal velocities of water particles and 
compared the results with the varieties of wave characteristics which have been studied by 
many engineers. 

With the use of the knowledge on shallow water wave theory, the horizontal component of 


particle velocity is expressible as 


U= A ain 2a (%, - Yr) bpd ert OPM MITC) 


N=("Vr)cosh (22 2/. )/sinh (21thfe,)s vom Se Anat) 


derived by the theories of Airy's small amplitude wave and. trochoidal wave. H,T,z, and L 


) 
| 
| 


where 


are wave height, period, distance from bottom, and wave length, respectively. 
For h/L<1/25, the waves are classified as long wave or solitary wave, and the maximum 


velocity of water particle is given by 


Usmex =aintH) * Yn sH) + piecing Balen 26) 


The wave characteristics, comparing the above theoretical values and experimental ones 
obtained in the experimental runs, will be treated. 


Fig. 6 shows the variations in maximum horizontal velocities of water particles at 
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pus points in Fig. 1 in the vertical 
tion, and z/h is taken in ordinate 
the ratio of water particle velocity 
at height z to bottom velocity upmax 
bscissa. The solid line in the figure 


ne theoretical curve computed from the 
y of small amplitude wave and trocho- 
wave as in Eq. (4). All the measured 
2s are the velocity of water particles 
rds shore. 


It is seen in this figure that at 


\ © @ Postion NOS 

rvation points where the relative fo) ooe aots » die odes 
© 2064 00218 oD © 00620 0.0218 
e depth, h/L, of incoming wave is in aD 0 ace2i 00378 @ 007398 aasze 


range of 1/lo~1/2 the experimental Solitary curve 


2s are rather agreed with results in 


shallow water wave theory curve, 

ite of various values in steepness. 
2ans that the motion of water particles Fig. 6 Relation between uyma /Upmax and 2/h 

be indicated by the theory of shallow water wave, even when the wave configuration of 

i waves characterized by small values of steepness in this range is smaller than that 
sng and solitary waves. At the observation point @), however, where the wave characteris- 
are in the range of h/L = 1/25~1/10, experimental results do not agree with the curve 
allow water wave theory, and the waves becomes similar to the long waves in behaviours 
the velocity profile in waves becomes uniform. Nevertheless, the experimental curve 
sates a poor correlationship with the curve derived by solitary wave theorys From these 
ances, it will be concluded the motion of water particle due to the waves in this region 
10t be treated by both theories expressed by Eqs. (2) and (4). 

Furthermore in shallower regions, no appreciable changes in water particle velocities 
»bserved in’ the direction of water depth. At the point @) however, where the water is 
shallow after passing breaking point, velocities near surface are greater than those 
bottom for waves with large steepness, whereas the velocity profile in waves with 

| steepness, is largely characterized by the uniform distribution. 

This evidence will result from the characteristics of breaking of waves. For large 
mess of waves, the breaking of waves is characterized by the formation of surface 

sr with white bubbles known as the spilling, and consequently the particle velocity 

“ds shore is larger in the upper layer than that in the lower layer. On the other hand, 
small steepness, the wave breaks suddenly as plunging breaker, it maintains the 
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characteristics of long waves and therefore the velocity change is not appreciable. 
Fig. 7 shows the relationship between the vertical distribution of maximum velocitigg 
of water particles taking z/hn in ordinate and the ratio of Uzmyxto Usmax obtained from Eq. 
in abscissa. This figure indicates the evidence noticed in Fig. 6 by slopes of empirical 
curves obtained by the experimental values. In the figure the solid line is the onshore ong 
In the foregoing discussion, the comparison of the theoretical values in small amplitude 
theory expressed by Eqs. (2) and (3) to experimental values is not referred, but it can be 
revealed that the experimental value exceeds the theoretical value at all points. The onsh 
velocity is closely approximated by the value in solitary wave theory in the range from shop 


side of transition region to near breaking point, as in the point (). This is verified by 


the results of field observation which Ke 


the Beach Erosion Board has performed at 08 


7) 
water depths of 5-15 ft. 


Fig. 8 is the relationship between 
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the dimensionless locations expressed by 
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h/Lo in abscissa and the ratio of maximum 


measured velocity of water particles at . 10420 eT = 
bottom, Upmax to theoretical particle — tae 
velocity by the solitary wave theory, 

Usmaxe The breaking line is also shown 

in the figure for reference. It is seen a A; | 


Position NO4 


that velocities at bottom approach grad- p04 
ually to those of long waves of finite a2 © © 00413 00103 
© ® 0052! 0026 
amplitude and the coincidence of both J nt ba Lassa oe oe 
10 20 
values in velocity is obtained at the ee 
‘San 


point a little near shore line from break- 
ing point. amax!Y omar 
(2) Variations in maximum velocities of 


water particles by locations. 


The fundamental ideas on sorting of 


= Ubmax/y), mox 


sand due to wave action and the theory of 
critical wave steepness on equilibrium 


slopes of beach are mainly based on the 


variations of maximum velocity of water 


particles along the beach profile and of 


mass transport velocity of waves. Using 


Fig, 8 Variation of upmax/v 
these ideas, a graphical representation Locations 


smax &t various 
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e relationship between the relative water depth h/Lo expressed in the abscissa and the 
| of maximum velocities of on- and offshore currents (uyzmex ANd Ugz7maxy Tespectively) to 
aximum bottom velocity upomx at a point of h = 50 cm in ordinate, (the maximum onshore 
m velocity and the maximum offshore bottom velocity are almost equal in their absolute 
s at a point of h= 50 cm.) is seen in Fig. 9. In the figure, the left one indicates 
elocity variations at bottom (z= 0), and the right one those at z= 2.2 cm. Generally, 
considered that larger sizes of sediment particles must be deposited near shore line 
the increase of difference between the maximum velocity of onshore current and the 

um velocity of offshore current. 


As seen in 
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endent of the 
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101 4 g mel hs 4 4 eee Le 
ring the 0 a0! 002 003. 004 Q05 206 003 004 005 006 

o/1 Oe, 

_ observation 
diment at Fig, 9 Variation of maximum water particle 


velocity at various locations 
al coasts that 


verage diameter of particles at bottom are the largest at points near breaking point 
he sizes of particle become smaller in the offshore. This fact will result from the 
tion in velocities shown in Fig. 9. For the waves with large steepness (H)/L,=0.0328), 
e region of shore side from breaking point, the difference between wzmax /Ugzmax and 
JYpomax at z = 2.2 cm is much greater than that at z = O cm, but for the waves with 
steepness (h/Ip = 0.0218, and 0.0103), the difference between Uy pmax! Ay maxed Mooney! Ubonex 
= 2.2 cm and that at z=0O cm are almost equal. 
From above descriptions and the distributions of mass transport velocity, the influence 
ich will be stated in (4) of this chapter, it is possible to clarify the difference in 
eteristics of current velocity with the changes in waves characteristics. The details 
be discussed in (4) of this chapter. 
Variations of water particle velocity in course of time. 

the previous sections (1) and (2), variations of maximum velocity of water particles 


e vertical direction at various locations are clarified and the sorting action of sand 
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drift can be discussed using 


the results obtained. As ha27” No2 


stated in the introduction, 
however, the change of par- 


ticle velocity in the course 


of time, namely the accel- 
eration of water particles, 
must be completely treated 
to establish the dynamic 
relationship of forces like 
drag and resitive forces to 
sand drifts. 

In Fig. 10, the time 


changes of water particle 


velocity at bottom of 


Velocity of offshore direction . Velocity of anshore direction 


observation points are in- 


dicated for the wave with 


d 


steepness of H)/Lo=0.0103, Fig. 10 Water particle velocity in the course of time . 
and the velocity changein the onshore direction is given in the upper domain above abscissa 
and that in the offshore direction below abscissa. Evidently seen in this figure, within | 
the range of shallow water waves, the onshore velocity possesses its maximum value at the 
crest of wave, whereas the offshore velocity its maximum value at the trough, forming a | 
sinusoidal variation. This curve is not completly sinusoidal and the slightly large aie 
on onshore side is indicated. The fact corresponds to the bottom velocity difference betwe 


on- and offshore currents in the wave theory of Stokes, which is indicated in Eq. (5). 


= an [RH Y coth (204 ){ coth [eah/p 1} 
Ati) Go ltr era re cen a oe ete, 


The ratio of the maximum bottom velocity of onshore current Umax to that of offshore 


current Us.,,. expressed by Eq. (6) is compared with the experimental results and illustrate 


eee 3 (mtLy ath (21/.) DOSS Sent 


roa Bids se, 
U2zmax, C a oe Sinh ath], * 2 Li sinh (gith/L) 


where C is the cerelity of wave. 


in Fig. ll. 


As seen in Fig. 10, the wave is transformed graduatly asymmetric in its shape when 
progressing to the shore line. The theoretical velocity given by Eqs. (2), (3) and (6) can 
be applied only to the shallow water waves in the range of h/L=1/10~1/2. Since the 
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‘celeration of onshore current increases as it moves 


{06 
om transition region to shore, whereas the accelera- 
O 00103 
on of offshore current becomes comparatively smaller, g ®@ 0.0218 
| eg (104-0 0.0328 
, is impossible to analyze the motion of water particle SS ~ 
this region, under the assumption that the particle Ss 
~~ 

tion is sinusoidal. Q 102 
| € 

The depth-wise distribution of acceleration is 5) 
irticularly considered to be one of significant factors 10 
1at indicate the movement of sand drift. Therefore, 1.0 1.02 1.04 1.06 


Experimental Value 


ie further theoretical and experimental studies to 


Fig, 11 Comparison of ohserved data 


stablish this relationship are distinctively needed. with computed results on ujmax/lemax 


+) Mass transport velocities. 

The variations of water particle velocities described in section (1), (2) and (3) have 
1 important influence on the instant motion of sand particles, and especially it is easily 
wWderstood that the water mass movement in one period will take an important part for analysis 
* net sand transport in the same duration. Eagleson, Dean, and Paralta’ have derived the 
slation between net velocity of sand movement and water mass transport velocity. 

As the first step for this purpose, the distribution of mass transport velocities in 
ater of 50 cm uniform depth in which the waves are characterized as shallow water wave, will 
> treated by the use of experimental data and theoretical results of Hamada, Sato, and 
Oeiect-Higcins. The theory of Hamada is derived by an extension of Stokes! theory in shal- 
ww water and the resulting equation is indicated in Eq. (7). Eq. (8) is derived by Sato, and 
1. (9) is the solution of boundary layer equation for velocity of mass transport obtained 


der the condition that the surface undulation is expressible as Eq. (10). 


Hamada's Solution: Giay=k? (LY: C: (cos kh /sinh kh): Sinhhz-cosh kB; 36 oo oe 


Sato's Solution: SS GREE Attar ahaa tary Uso) 


Ponecuy Meetnees ee - s akh-f 3(Z Jerk 
Siena U(z= z h 
Solution: fe aa 3 (Sanhzkb o Pa steals Ae EE ee 
ey t my | ) Weel 
Sette 24 a aparece, Mathys GO) 


ere (12) 


t+T 
GQ= ae oa 2/7 Gia-(/7): J, udt Si) ace Hoa 


vertical distributions of mass transport velocities for the offshore waves with the 


ues of 0.0328 and 0.0103 in steepness at point @) are shown in Fig. 12 and it is seen 


at Eq. (7) indicates a good agreement with the experimental values for the wave of high 
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steepness (0.0328), whereas Eq. 


(8) is a better expression for the 


14 - Q.0/03 
' 


wave of low steepness (0.0103). 


Consequently, it will be concluded 
that for steep waves (Hp/Lo= 0.0328) 


the mass transport velocity ap- 


proaches zero at bottom and rapidly 
increases in the upward direction, 


whereas for gentle waves (Hy/Lo = 


0.0103) the onshore velocity is -05 O 05 asc -O5 Oo O5 I pac. 
appreciable in a form of nearly Offshore dredion Onshore direction Offshore direction Onshore chrection 
uniform distibution. Therefore, Fig, 12 Variation of mass transport velocity in 

vertical direction 
Longuett-Higgins's equation derived 
under the assumption that the total mass transport is zero at all cross sections and the 
wave shape is represented as clapotis as seen in Eq. (10) can not agree with the experiment 
data obtained in the region where the backwash is not predominant. 

The depth-wise variations in mass transport velocities at onshore side from the transi- 
tion region are next concerned. Values of mass transport velocities at bottom in the region 
where it is pressumed that the influence of undertow after breaking becomes appreciable are | 
zero when calculating by Eq. (7), whereas positive in Eq. (8). These two equations, con- 
sequently, are not suitable for computations of net mass transport velocities including the 
mass transport by the undertow. 

Now, assuming that a constant and uniform back current is produced in the whole depth 
the velocity U, of back current is expressible by the continuity theorem as follows. 


L h 
(Yat) dx | (ut Us) 4z =0, oo ore eee Gee 


h 0 
Ge= (Mn) Uw) dz = (aR C/2h)- co kh, ee eo ee 


where U(z) is the mass transport 


velocity expressed by Eq. (8) and 


a=H/2. Fig. 13 shows the compar- 


ison of experimental values to 


computed values by Eq. (8), putting 
U(z) =U, at z=0 and U, - Up=Uo 


and Cy, the wave cerelity in deep 


water. It is seen that both curves 


; P Sig, 15 Variation of mass transport velocity at. 
are similar in their behaviours, bottom 
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therefore it is found that the velocity of undertow shown by Eq. (13) must be included 
omputation, in order to clarify the flow behaviours in onshore side of breaking point 
h is largely influenced by the undertow. 

The vertical distribution of mass transport velocities at a point, where the influence 
ndertow becomes appreciable, is expressed as follows, using Eq. (8) and Eq. (13): 

HY only, HH), coshate/L) H coth (2A(L)}-C . ears) 
Cc Sh Toth] he 

Fig. 14 indicates the computed values 

Eq. (14) for the vertical distribution 
ass transport velocity with the experi- 
al data, and it is evidently seen that 
od agreement is obtained. 

As the results, the onshore current 
eds the offshore one at the surface of 
r, while the latter exceeds the former 
he bottom, in the shore side of break- 
point, for all incoming waves. This 


ence has already noticed in the report 


field observation of Eaton, Inman, and 


) EE Oe ee MR 
3 
7 Ux/e, x 10 Uk/c,x10° 


Fig, 14 Vertical distribution of Tt{z)/€, 


‘The influence of wave characteristics on U'(z) by the use of Eq. (14) is indicated in 
15, where values for U'(z) are plotted in solid line for gentle wave of H,/L,= 0.0103 
in dotted line for steep wave of H,/L,= 0.0328, when h and H, are constant. It is seen 
ais figure that the difference between U'(z) at surface and that at bottom is not so 
inct for gentle wave, while for steep wave, there exist remarkable onshore currents at 
ace and offshore currents at bottom. The states of currents, therefore,are mainly in- 
nced by the steepness of breakers, and it will be the basis of Hayami's idea on critical 
Be scpnasa. The water depth z, at which U'(z) transforms from positive to negative 


> it computed by Eq. (15): 
cosh 2R2.=(/2kh)-Sinh 2kh, i Oe Gh) 


Consequently, z/n is little influenced by h/L and 0.56~0.60 for all values of H)/lgas 


‘in Fig. 15. 
‘In the present chapter, the behaviours of water particle in wave motion, being mainly 


ential on the movement of sand, have been discussed. The movement of sand particles on 
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beach bottom will be treated by means of the preceding 


concept of mass transport velocity in wave motion. The 


meme 
aE =QQ/03 


° 


experimental investigation described in the foregoing, 


He 
“L, 70.0328 however, is still preliminary, and the complete clar- 
ification of the hydraulic behaviour in wave motion is 


not made. After the real success of the present exper= 


imental study is obtained, the dynamic relationship of 
—/ sand movements in terms of the introduction of mass 
transport velocity will be certainly established. 
e cts 0 o on “(2 However e engineering need arising from e 
Fig, 15 Effects of Ho/lo on U (z) in » th ineeri d arising from th 
case of constant Hy and h 
rapid progress of coastal works urges also to make a 
relationship of sand movement in dimensional or empirical fashions, so that the transporta= 
tion of sand drift will be concerned in the following, as an attempt being derived under 


the assumption that the movement of sand drifts in wave motion will be similar as that in 


’ 
open channel flows. 


4. Remarks on mechanisms of sand drift motion 
The evaluation of ue in waves which are mainly characterized by the acceleration term 
has not been made at all. Hence, the emperical relationships between the critical frictional 


na 


velocity ug* and JsgD in open channel flows will be used for the reciprocal movement of sa 


drifts in wave motion as engineering approximation. Referring to the study on the critical 
) 

condition of sand drifts by Shinohara, Tsubaki and Yoshitaka, the mean value of ratio of 

to JsgD is approximately 0.235, and the ratio of representative flow velocity acting on sands 


to the critical frictional velocity is 8.5. Therefore the following equation is obtained; — 


Uc/VAeD = 185, se ee nen 


where s: specific weight of sand in water, g: acceleration of gravity, D: diameter of 
sand particles. 

Table 1 describes the values of critical flow 
velocity for various sizes of sand grains. 

When the water particle velocity caused by the 
wave action exceeds the critical velocity, the sand 
initiates to move. In the behaviour of sand movement, 


of essential significance is that the relationship 


between the amount of sand drift and the motion of 
water particle must be established. However, no Table 1 Critical velocities for 


articular sizes of sand grains 
studies on the amount sediment transport in unsteady Z oi 
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ave been made and futhermore 


x07 


>rical formulas also have been 
1 in explicit forms. Therefore, 


Imption is made here that the 


-al relation for sediment 


ort in rivers applicable to 


»w in wave motions, and thus 


[lowing equation is used for 


tion of the amount of sand 


2} 2 
alto <A(u UE) + +» 7) 


is the amount of sediment 7 OMAN 
1B Fig, 16 Relation of u(u® - u,*)/C,5, h/L, 
srt per unit tine, u* is and o in case of D = 0,5 mm 
tional to u as stated before, 
is a numerical constant and approximately unity by experimental data, when (ux? = u,*) 


ymed not to be large!” Eq. (17) can therefore be rewritten as follows, for a partic- 


ize of particles. 
Z iz 
ge Gul(u-Uc), ese 


snoting the average amount of drifting sand in one period by ap> it can be expressed 
t+T 
3- (Yr) Sgdt oc U(uzUe), roean et aodeGeg) 
1@ average amount of sand drift at a point can be estimated, if u(u2 = u,?) are obtained 
Bis locations along the beach profile as shown in Figs. 16 and 17, which indicate the 


mental results for D= 0.3 mm and 0.5 mm, respectively. 


ue experimental wave steepnesses 


*e 0.0328, 0.0218 and 0.0103. 


% 16 and 17, the positive slope 


ses the process of deposition, 


3 
ry 


she negative one that of erosion, = 

; the sand ridge will develop at a 

; where the sign of slope 3 

rms from positive to negative. 

ie discrimination between depos- 

nd erosion regions in the beach 

-is therefore made by the use of Fig, 17 Relation between u(ue = u,2)/c,5 h/g 


5 OL 2) = O25 
6 and 17, and it is indicated in and Ho/L, in case © > 
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Fig. 18 for the case of D=0.5 mm. As a 
x 


——— 


D=05mm 


' te 


Region of Erosion 


seen in Figs. 16 and 17, the sand ridge 20 


is developed only by the steep wave of 


H,/Ly= 0.0328 and its location is in 10 


the region of h/L,= 0.03~0.05. On the 
other hand, when the incoming waves are "» 
gentle (H)/Lo= 0.0218 or 0.0103), the IN“ S 


= 
{ 
sand ridge can not develop, despite of S 
Ss 


all particle materials in size, and fur- 10 


| 


thermore the sand moves towards the \ Region BREOSiOn 


shore line, which may explain the grad- 
ual progression of shore line. The 


location of sand ridge formed and the 


beach deformation estimated by the 
dynamic relationship of particle move- 


ment in waves can be verified by the 0 
0020 


experimental results of equilibrium 


Aly, 


beach profile obtained before, as seen Fir. 18 Beach progess estimated by the variation of 
P , u(ue = ug 8 at various locations in case 


aboy lakes alleys of D = 0,5 mn, 


Adding to the above study, the 


[BOGE 


estimation of beach process in deposi- 


tion and erosion is described in Fig. 
20, under the assumption that the amount 


of sand drift is proportional to the ef- 


fective water particle velocity (u - ug), 


in which 


tT 
(a= (Yr) (U-Ue )dt + + (20) 


In Figs. 18 and 26, it is seen all Fig, 19 Equilibrium beach slope for D = 0,5 mm 


curves for beach process are similar in 
both cases. The available formula for estimation of the amount of sand drift is not obta: 
up to the present time. However, the mechanisms of beach deformation only studied in conn 
tion with the steepness for deep water waves can be sufficiently explained by means of thé 
idea of water particie velocity subtracted by the critical velocity Ug which is mainly in 


fluential on the movement of sand particles. 
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1clusions 
1e velocities of water particles in off- and onshore regions, surf zone and shore side 


uking point were measured with the use of the current meter of unbonded gauge type to 


igate the mass transport velocity of 


The actual behaviour of undertow in 


ore side of breaking point was also Region of Erosion 


ied. Although the mass transport 
cy is mainly influential on the 


it of sand drifts, the present study 


‘eliminary project of the whole pro- 
1 the dynamic establishment of move- 


> 


" sand drift can not reveal the basic 


Erosion 


> prénciple, so that the authors 
1 attempt to estimate the beach 


ation in similar forms of bed load 


ons in alluvial channels, with results Fig, 20 Beach process estimated by the 
variation of (u-ug )/Cg at various 
larification in beach deformation locations in case of D = 0,5 mm 


8 


ye authors intend to obtain the comprehensive solutions for mechanisms of shore line 

1 by analyzing further the complicated behaviours of time variation of water particle 
cies in the onshore side of transition region where the sand drift movement becomes 

2, of the critical flow velocity of sand particle in fluid characterized by continuous 
sreciable change in accelerations, and of the amount of sand drift transported. 

> is a pleasure to record here the authors' gratitude to Assist. Professors Y. Iwagaki, 
thara, and Y. Iwasa as well as Mr. J. Yamaguchi, for their very useful suggestions, 


isms, and cooperations. 
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Introduction 
In 1957 and 1958 the present authors investigated the suspended littoral drift at 
nakomai, where the industrial harbor is now under construction. The city is located on 

> southern coast of Hokkaido as shown in Fig. 1. The suspended sediment in the nearshore 
2a was captured by bamboo samplers which have been employed 
the authors in past several years for the investigation of 
)2) 


1 
11d movement along the coast of Hokkaido. 


A certain number of bamboo poles were vertically set in 


2 sea for a few days. They were drawn up carefully, so 


at the sediments captured in the bamboo samplers would not 


ow out from holes. Samples of sediments were then packed 
polyethylene sacks. The sediment was desiccated, and 


ighed in the laboratory. 


The time intervals of setting were as follows: 


19573 from 10.00~12.00 A.M. Aug. 2 to 2.10~12.00 A.M. 


Aug. 7 (120 hours) 
19583 from 10.00~11.30 A.M. Aug. 27 to 12.00 A.M.~2.00 P.M. 
ie eee al 
Aug. 29 (50 hours) 


2 distribution of samplers are shown in Fig. 2. The mean tidal range in this area is 


230 
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13 
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about 1 meter. During the each interval of observations, the wave heights were respective- 


ly 0.8 m~0.4 m (1957) and 1.0 m~0.3 m (1958). 


2. Vertical distribution of suspended sediment. 

As yet few investigations have been made about the vertical and lateral distributions 
of the suspended sediment in the nearshore see Some examples of the vertical distribu- 
tions of the mass of the sediment load captured by samplers are shown in Fig. 3. The 
relation between the mass and the height of hole from the sea bed is plotted on semi-log 


sheets. 


In most cases each curve has a bending point which shows that there exists a 'Sprungs- 
schicht' of suspended sediment about 1 m above the sea bed. In the upper layer the increasell| 
of the mass per 1 m (height) is much less than in the lower layer. The absolute value of 
the mass is also larger in the lower layer. Some curves, however, show that the suspended 
sediment is somewhat decreased near the surface. Perhaps on account of tidal range and 
wave motion of the surface water, holes at the upper part of the sampler can not always 
capture the suspended material. Furthermore, they are sometimes exposed to the air and 
lose the material which is once captured by bamboo poles. 

Grain-size analysis was made by sieving the samples in order to study the size distrib-— 
ution of the captured sediment. As will be seen in Fig. 3, the suspended sediment in each 
depth contains various grades of size in a similar proportion. This fact is nothing but 
a feasible support of the actual existence of the double layer. 

The vertical distribution of mass thus obtained are classified into the following 
four types which are shown in Fig. 4. 

Type I. This type of distribution is understood to be the typical one which appears 


in natural sandy coasts. The double m9 


O>.40mm. ©>.1 1 mm. 
©>.25 8>.075 
©>.149 @<.075 


layer, mentioned above, can be seen 8 


! 
¢ 


in this type of distribution. The 
slopes of the curve are 0.23 and 7.60 


per meter of the depth in the upper 


and lower layer, respectively. Some- 


times the mass decreased suddenly at 


the surface in the upper layer. 


Type II. This type appears at 


es ! 
SN | TT RS 
SUITE TTI siete | TUT 


.Ol 100 
ergy of waves concentrates here, the gr 


the corner where a’ sandy coast and a \ 


breakwater meet together. As the en- 06 
turbulence of sea water is quite Big. a5 
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irkable and the curves of the 


m 
sar layers are steeper than in a all 
8 
rot Hill 
sr layers is somewhat convexed a # ASEH 
ry) 
the right showing that the sedi- ea 
t load is poured into the corner a5 
ul 
ugh the samll river indicated in 4 | 
2 il 
> ° 03 t) 
Distribution curves in this 4 hook 4 iil 
me 5 0 
eeeemeanipeeemeeg | MRIIASUIL; UI 
in | 4 
2d with type I. If an observation = MB ilitestaa Ea ~J 
i 


(e) 


made at several points between 100 ar. 
sampler of No. 1 and the shore- "Mass OF SAND DRIFT CAPTURED IN BAMBOO PIRES 

2, the curves in this zone must be Fig. 4 

ilar to the distribtuion of type II. 

Type III. Five curves observed at points of No. 16~~No. 20 belong to the distribution 
sype III. They do not have the lower concentrated layer. These points are located at 
corner formed by breakwater as shown in Fig. 2. As compared with the other types of 
same depth, it may be understood that the concentrated load in the lower layer of the 
s III is protected from the disturbance in the more shallow zone along the shoreline. 

Type IV. The absolute amount of mass captured is generally less than in the other 
2s. The upper part is somewhat vertical. The development of the lower layer can be 
2 as in type I. 

Points of the observation where these curves were obtained, are located at the west 
: of the harbor. The coast in this area is now being eroded gradually. 

According to G. E. Vincent it is reported that, an entrainment current which transports 
sediment load develops along the direction of mereae Shinohara and others noticed in 
ir experiment of suspended pulverized coal that a thin layer was moved by waves with a 
sp slope of vertical concentration in the immediate vicinity of the pea” They called 
layer as ‘transition zone'. In the case of sand, this layer was not observed. In the 
d observation which was carried out by the present authors, the velocity and the direc- 
1 of the bottom current could not be determined. Therefore, it is still questionable 


she lower layer of heavy concentration represents the entrainment current proposed by 


ent. 


Lateral distribution and profile. 


The integration of the suspended load from the bed to the surface were obtained for 
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each point. The results in Fig. 5 shows the lateral distribution of the logarithm of 
these total mass captured in a day. The total amount of the sediment at each point de- 
creases geometrically with the distance from the shoreline. The total amount of the mass 
in one point depends almost upon the mass in the lower layer excepting points No. 16 ~ 
No. 19. In this meaning, Fig. 5 represents the degree of development of the lower layer. 
By using the curve of type I, the equi-concentration curves of suspended load along the 
profile, i.e. normally to the shoreline, are obtained in Fig. 6. The depth was observed 
by lead sounding. Profile obtained in 1957 clearly shows the relation between the sand 
ridge and the concentration. Inside the ridge, the concentration increases rapidly to the 


shore. The ridge was not remarkable in the profile observed in 1958, because points of 
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observation were selected much nearer to the harbor and the number of points was only 


half of the ones in 1957. In 1958 the observation of the lower layer was not successful. 


The authors are very much indebted to the stuff members of the Laboratory of Applied 


Physics, Faculty of Engineering, Hokkaido University. 
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1. Introduction 

Primary considerations of a harbor engineer in planning harbors on gandy coasts are 
usually the following: 

a. A well planned harbor on a sandy coast should require little or no maintenance. 
Therefore, the first consideration is how to plan a harbor relatively free from excessive 
deposits of silt (accretion) thereby resulting in a low maintenance cost of the harbor. 

b. A harbor well protected from drifting sands is apt to create a problem for ships 
entering and leaving the harbor. Hence, a compromise will have to be made to satisfy the 
two conflicting problems. 

c. The bigger the harbor is, the more it becomes subject to accretion and erosion 
of the adjacent coast-line. It is often the case that the land value of such coasts are 
often high in comparison to the scale of the harbor. Therefore, consideration must be 
given to the use of accreted lands and the counter-action of erosion in the planning of 
the harbor facilities. The size of the harbor most suitable to meet these conditions 
must be considered. 

The problems confronting a harbor engineer are conflicting in nature, but must be 
solved in a rational manner. In order to solve these requirements, the following data 
should be looked into for present-day coastal engineering: 

(1) Predominant direction and quantity of littoral drift. 

(2) Dimension and frequency of prevailing waves. 

(3) Relation between water depth and drift amount. 

(4) Optimum side slopes and the advancing and retreating of the shoreline. 

(5) Design and engineering data of adjacent and other harbors. 

Item (1) is sufficient to solve item (a), i.e. to decide the location and direction 


of breakwaters and port entrance, and then to work out the amount necessary to be dredged 


within the port. 
Item (2) is a fundamental consideration in planning navigable channels, and for esti- 


mating wave movements within the harbor. Item (1) of course must be considered for selec- 


ting the disposal areas of the dredged material. 


-—-59- 


THE MEASUREMENT OF LITTORAL DRIFT AT SHORELINE HARBORS. 
\ 


Item (3) must be considered in order to determine where maximum drift will occur 
along the ccast-line, thereby permitting selection of the harbor entrance so that there 
will be a minimum of drift into the harbor and shoaling of channels outside the port. 

Item (4) must be considered in order to estimate the scale of accretion on the front, 
and amount of erosion at the rear of drifting sand around a breakwater during certain 
economic years. This consideration must be made in addition to Items (1), (2), and (3). 


Item (5) will be a compilation of data and past experience. 


2. Geographical survey 

This section is devoted to study of deriving at a method of judging the phenomena of 
drifting from geological and topographical viewpoint. Needless to say there is a tremen- 
dous amount of energy dissipated on coast lines, and the yearly amount of wave energy spent 
on a simple coast line facing the open sea for a coastline of 1 km is 107' erg, comparable 
to the energy of typhoon which is on the order to” to 10” erg. Littoral drifts have gross 
movement in addition to the irregular effects of the winds and waves, and it is not an 
easy task to make analysis from short-term observations of data. 

Common methods currently adopted in determining experimentally the direction of drift- 
ing sand from geology and topography are the survey of transformation of shore patterns 
in the vicinity of heedlands and inlets, the measrurment of sand size distribution of bed 
materials or classification of ores, etc. 

The author has made a study of the Yufutsu Plain (the backyard of Tomakomai harbor), 
from the time of its formation from the point of geomorphology and then assumed the condi- 
tion of the littoral drift around this area. 

The island of Hokkaido is considered to consist of two parts; the Ishikari-Chitose- 
Yufutsu line was a strait having changed into land during a geological year; this land is 
called the Yufutsu and Ishikari plains. A hill zone north of the Yufutsu Plain was formed 
from a land-rise, but the birth of Yufutsu Plain itself, involving vast marsh lands and 
composed of sand only is considered the result of littoral drift. 

An aero-photograph shows big are lines of old sand-hills, and that is the chief basis 
for assuming the above process. Furthermore, earthen-wares and shells dug out from shell 
mounds around the Yufutsu Plain serve to prove that the back arch line of sandhills, that 
is, the old coastline is of about 4,000 years old. (Different kinds of shells are known 
to live according to the amount of salt in the sea. Therefore, it is easy to determine 
from the location of where the shell was dug out whether it was sea-bottom or land before.) 
It is concluded that 4,000 years ago littoral drifts caused by the Hidaka River deposited 
in the so-called "Yufutsu Bay", shoaled and made narrow there, advancing the shore-line 


with some parts of bogs left, and then about 300 years ago they formed same shore-line as 
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Photo-1 
le present one. For this reason we can find a layer of volcanic ashes caused by eruptions 


#7 and 289 years ago. (The age of a layer of volcanic ashes may be exactly presumed by 
sasuring radioactivity of an isotope C-14+ from grasses in the layer.) In this way, the 
ithor has tried to develop his opinion of drifting sand under the following assumptions: 

(i) The accretion of littoral drift around the Mukawa-Noboribetsu low surf zone, in- 
uding the Yufutsu Plain, began 4,000 years ago. 

(ii) The amount of littoral drift has been constant for 4,000 years. 

(iii) There has been no change in sea level. 

(iv) The present shoreline was formed about 300 years ago, and no change has taken 
ace after that time. 

(v) For this period all drifts along the shore less than 5 meters in depth were im- 
unded and accreted. 

According to the above conditions, about 45,000 cubic meters per year of drifting 
nds along the Tomakomai coast are estimated. This figure is considered to be fairly ac- 
rate, judging from the amount of 14,000 cubic meters of annual mean sediment at the up- 
ast side of Tomakomai Harbor, which has been calculated from sounding survey maps for }+ 
ars. The prevailing direction of drifting sands along the shoreline is from East to 


3) 
st, since estuaries in this area are all facing West. 


Measurement of forces moving sand 
) Wind: It is not easy to make long-period observations of the waves having the closest 


nnection with drifting sands. It is, therefore, convenient to study the wind which is 
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most closely related to the wave, and relatively easy to make recording of such wind 


readings. Winds, on the other hand, are greatly influenced by nature and due considera- 


tion must be allowed. 

Let us now look into the case of Tomakomai Harbor. There are 4 observatories called 
Hachinohe, Muroran, Tomakomai and Urakawa, surrounding the off-sea side of the harbor. 
First, the wind data at the above observatories were obtained, and then a wind velocity 
frequency diagram, a wind direction frequency diagram and a velocity-direction distribu- 
tion diagrams for more than 5 m/sec, 10 m/sec and 15 m/sec were figured respectively for 
one year. Next, a study of wind characteristics for each season was made from the 28-year 
data at each observatory. With these characteristics, a trial was made to see how strong 
winds influenced the tendency topographically. As a result of this study, in judging the 
wind off the Tomakomai coast, the record of Urakawa Observatory will have to be used for 
Westerly wind directions, Hachinohe for North-Esterly, Muroran for South-Westerly, and 
Tomakomai for South-South-Easterly winds. The results coincide with the forecast made from 


the point of topography. 
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Next we ascertained that the 10-year records the Tomakomai Observatory had no re- 
kable change for a period of 1 year and then worked out the maximum wind velocity, and 
maximum duration time in each direction for 2 years. Thus we arrived at the conclu- 

n that all winds without S-SE direction were not relied upon, and that the wind off 
akomai Harbor was predominantly from the East in the summer. In this way we could 
© presume the prevailing direction of drifting base 
The effect of waves: In order to solve sedimentary movement, it is necessary to 
erstand statistically sea roughness, instead of accurate observation of each wave. 
refore, we will have to deal statistically with-the results of a short -time continu- 
observation, or a regular intermittent survey, and then have to summarize the results 
© a whole-year sea roughness. 
Being engaged in the planning of Tomakomai Harbor which is situated on a sandy coast, 
- author found out that a wave recorder operating at a depth of 9 meters one year, from 
e 1, 1955 to May 31, 1956, gave the statistical distribution of wave height and period 
the form of log-normal Cats & Proueel lity diagrams of wave height and wave period are 
wn in Fige2 ana Fig-3 respectively. 
The equation of probability density of each wave height and wave period are as follows: 
(a) Wave Height: 
Y= \/(le7 x 2.2844)exp - (bog -/. 760) (2x 02544) 
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(b) Wave Period: 
(1) For fetch infinite direction 
2 
Y= \/GeR x 2425) 2xp -(X-8, 627)/( 2% 2.425) 
(Hi) For fetch definite direction 


y = \/(2% - 0.005 ) an -(X-4/69)/(2x 005) 


Probability density 


where 


4 


x 


Significant wave height and wave period in each equation 


In this way, the probable percentage for waves of more than 4 meters high plotted 
against time is 0.15 %, for 4+ to 1.2 meters - 13 %, and for 1.2 to 0.4 meters - 58%. As 
for wave period, there can be divided into two groups of 4.17 sec and 8.63 sec and the 
occurrence probability complies with the law of log-normal distribution, being a function 
of driving distance. The same can be said for wave energy occurrence probability. For 
the present time observations have been made regularly twice a day for period of one year. 
In — to study littoral drift, not only should further precise values of various sta- 
tistics be made, but also more frequent measurements be made, i.e. 3 or 6 times regularly 
daily for a few years period, are needed. | 

As shown from the above, the effect of waves on a shoreline, and furthermore by cal- 
culating the wave energy in each direction and summarizing as shown in Fig.4, the prevail- 


ing direction and the scale of littoral drift can be predicted. This information has 
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1 valuable in determining the normal line for planning the breakwaters against the ef- 
“s of sand and waves. 

The effects of currents: In general currents along a coastline can be divided into 
ul current, ocean current, density current and wave flow. The first three currents are 
general constant and are not of a great magnitude for each given area, with a few ex- 
tions. Currents caused by waves have a tendency to alter its direction and velocity 
sanding on the condition of the wave. The current having the greatest effect upon lit- 
al drift is this one caused by waves, and is generally called “littoral current". Hence 
study of currents is very complicated, and in many cases records obtained after pains- 
ing effort must be discarded as unreliable. Though many formulas have been published 
calculating coastal currents based on winds and waves, none seems to be adoptable. 
le 1. shows a wave-current relationship at Tomakomai area. Since the coastline at 
akomai runs in a East-to West girection and faces South, littoral currents caused by 
2s from ESE - SSE and S - WSW are driven from West and East respectively. According to 
table, littoral currents in a stormy weather, which has the greatest effect on sand 
“ting are entirely controlled by waves; and the bigger the angle between a wave direc- 
1 and the shoreline, and also higher the wave, the stronger the littoral current. 

Measurement of slopes of sea bed: In order to predict the process of accretion or 
sion caused by breakwaters, it is of importance to know the slope in equilibrium of the 
st line. Roughly speaking, a bed slope can be given as a function of wave elements and 


in diameters of the bed sands. Consequently, at coasts where neither waves nor bed 
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sands change from year to year, this can be performed by measuring the bed slope with 
extreme care. Especially, it is important to take measurements of accretion and erosion 
around structures (under construction) than it is on a natural coastline, and taking an 
Keardce on a yearly basis. 

Fig.5 is an example of changes of slopes of an accreted area on the front side of a 
breakwater under construction at Tomakomai Harbor. The length of the breakwater was ap- 
proximately 60 meters, and the foreshore slope above low-water-level was almost constant. 
There was no effect of accretion at distance of 70 meters toward the seaward side of the 
breakwater. The same tendency can be seen at Irako Harbor (Aichi Prefecture) which is 


also located on a littoral coast. 


5. Measurement of littoral drift with isotopic tracers: 

There are many methods of measuring littoral drift, and one method tried by the suthor 
7), 8, 7), 10), 11,12), 14). 
with the use of radio-isotopes is described herein. ¢ 


This method consists principally of the following; prepare homogeneous radioactive 


glass beads having the same specific gravity, and the same size distribution as that of the 
sea bed sand, and measuring the effects of coastal currents, waves, etc. with radioactive 
tracers after placing the glass beads on the sea bottom. Actual results of these measure- 
ments made at Tomakomai Harbor by the author were published in the "Dock and Harbor Author-_ 
ey! Mand will not be repeated here. Similar measurements have been made in England at the 
Thames Estuary, Poole Harbor on Norfolk coast, and also at Kaike Coast Tottori Prefecture, 
Irako Harbor Aichi Prefecture in Japan. 

The radio-isotopes used in Japan were Zn-65 (Kaike Coast and Tomakomai Harbor), and 
Co-60 (Tomakomai and Irako Harbors), that used in England was Sc-46. Sc-46 excells in 
gamma-ray flux density and sensitivity to counters, but has the disadvantage of shorter 


life, being 85 days. If an atomic pile is available, Sc-46 would be the best medium to 
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It would be advisable to make the diameter of the glass beads uniform and approxi- 

the majority of the diameter of the sand granules of the bed sand, thereby simplify- 
analysis of the test results. The grain size distribution of the radioactive glass 

S was made the same as the bed sand for the experiment at Poole Harbor, as in the case 
omakomai, but according to regulations of the Radiation Hazard Division ot the Harwell 
ratory, those beads with over 0.1 micro-curies per grain (or over 0.25 mm dia.) were 
ved. The glass beads used at Tomakomai were less than 0.2 micro-curies, and so there 
no radiation harzards. Care must be taken to place the glass beads on the sea bed in 
ng with the bed sand. In the experiment at Poole in England, the glass beads were drop- 
from a height of 8 ft (depth of water 19 ft) and it was found that the glass beads were 
ersed by the ocean current and making a movement different from the bed sand, that from 
second time they made it a practice to release the glass beads on the sea bottom. At 

e glass beads in the amount of 1.75 curies per observatiom were used and it is claimed 
they were able to make continuous readings for a period of 4 months even during the 
hest winter months. On Norfolk Coast in England it is claimed that they were able to 
continuous readings for 3 successive months at a depth of 5 meters. At Tomakomai it 
possible to take readings for a month and a half for a amount as small as 2 to 30 milli- 
ss. It seems that the radiation counters of Japan are just as effective as those of 
and. It goes without saying that the efficiency of the radiation counters must be 
ed before use. This includes the directional ability, dead-time, absorption and back- 
tering of the sea water and bed sand of the detecter. Fig.6 gives the radiation ab- 
tion of the gamma-counter for cosmic ray measurements, which shows that the count of 
10,000 counts falling off to 10 counts at a distance of 70 cm. It can be deduced 

this that the radiation count is limited to 
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covering up of the glass beads. The resluts of the experiment at Tomakomai can be summed q 
up as follows: 7 
(a) At depths of 3 to 4 meters, and near the shoreline, the movement of bed sand was 
so great as to not be able to give relationship of the wave to drift. Improvements must 
be made in the method of measurement. : 
(b) At depths of 6 meters, no perceptible movement of bed sand is noticeable against q 
waves of 1 meter in height, and comparatively small granules move for waves of 1 to 2 
meters in height, and largest of sand granules move against waves over 2.5 meter in height. — 
The sand and waves move in the same direction, with the amount of bed sand moved in propor=—_ 
q 
; 


tion to that of the wave energy imposed upon the sand. 


(c) At depths of 9 meters, bed sand movement is not noticed for wave heights of 2 


meters and under. Any movement of bed sand is subject to the ocean current. The bed sand ~ 
q 
showed some movement against waves over 2.5 meters in height, but the direction of the move- 


ment did not agree with the wave, and the direction of the ocean current had a marked influ- 


ence on the direction. j 
(a) Bed sands near underwater structures (breakwaters) are controlled by waves rather 
than by littoral currents caused by the structures themselves. ; 
(e) In summary, it can be said that the movement of bed sands is more influenced by ) 
waves than by littoral currents. ; 
Measurements have been underway since January 1958, by radio-isotopes at Irako Harbor 
in order to study the shoaling process caused by littoral drift. To date the following ‘ 
has been eee : 
7 


The grain size of sand and the regional conditions of external forces has allowed us 
to determine whether the sand movement is of the creeping type or the suspending type. 
Fig. 7 gives the two modes of dispersion of radioactive sand. The contour lines in the 
figure indicate zones in which equal numbers of radioactive impulses were received. The 


resulting diagram thus gives the dispersal of active sand around the initial injection 
point. 


6. Items of measurement of littoral drift and harbor plannings 

(1) Items to investigate: In making preliminary studies of harbor planning, the 
first consideration is to make a geological survey, and then make a rough determination 
of the items given in the introduction 1. and consider these items in the light of observa= 
tions of the phenomena taking place on the coastline at present. : 

In making studies of the wind, it is of prime importance to obtain wind data from 
Observatories located in spots where they can accurately catch wind date creating waves 
affecting the proposed harbor construction site. Many cases are cited where due to insuf- : 
ficient wind data littoral direction has been misjudged. We believe that next to geologi cll 
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| R.1. SAND INJECTED MOVED BY WAVES, ETC. 


CREEPING TYPE 


R.1. SAND INJECTED MOVED BY WAVES, ETC. 


isiderations, wind study should be placed topmost in making coastline studies. 

After preliminary studies are made, the next step is to make full scale studies. At 
i's stage the most important factor is the study of waves. Since littoral drift is in- 
ienced mostly by wave energy, it is of utmost importance to know the elements of the 
res and the density of occurence. Though automatic wave recorders are still to be 
‘fected, the only way to study waves is survey waves from high places, or by making a 
‘raction diagram each time. It seems that a study of currents is a waste effort as far 
littoral drift is concerned. In general, the changes in ocean currents are as large as 

the coastal currents and tidal currents for prospective harbor sites. Hence, a one=- 
‘e observation will not yield satisfactory results. Since currents are combined of wave 
rements (mainly ocean currents) and irregular coastal currents, the changes in currents 
» unstable as the winds. Since satisfactory recordings cannot be made of currents this 
dy is not recommended too highly in making harbor studies. 

The measurement of bed materials is significant in relation to bed slopes in equilib- 
im, especially predicting the slope due to accretion and erosion after completion of 
akwaters. It is important to measure the slope of accreted deposits while the construc- 
mn of underwater structures is underway. This study will not have to bid during rough 
ther and if geological surveys have been conducted, the slopes during rough weather can 
determined. Since it is not easy to study the bed bottom during rough weather, such 
dies are not necessary beyond the immediate vicinity of the structures. 

Measuring with isotopic tracers has only begun recently and has given engineers a big 
erest because they can make a direct measurement of drifts for a long term. Direction, 


tance, or velocity, except for amount of radioactive sand will be clearly recorded even 
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if some were to be covered up. As to the amount moved, we have to expect for a pattern 
experiment. Anyway various records of measurement of Saceeat phenomena will not be use- 
ful unless they are arranged at least to show annual tendency. A statistical arrangement 
is a key point. 

(2) Harbor planning: Any good survey of a natural coast line will not give an exact 
knowledge about the processes of accretion, erosion and shoaling after building-up of a 
harbor. 1t is, therefore, necessary to look into examples of harbor construction in the 
world and to produce a key plane Then construction can be started in a "try and check" 
way rather than in an "all-over" way. 

Drifting sand will come in a harbor on a sandy coast, however it may be built with 
reasonable planning. But if its entrance dees not face the prevailing way of drifts, it 
will be pretty free from shoaling. On the other hand, the site of the entrance is mostly 
decided by the site of wharfs or the scale of ships coming in and out rather than by 
considering the depth where there is little drifting sand. 

In deciding the position and direction of the entrance, it is of value to estimate 
what ithe filling pattern of deposits might be expected at a proposed nearly barrier. 

Fig.8 is an example simply sketched of Tomakomai Harbor. The assumptions are: The edge 
line AB of the deposit area is perpendicular to the direction of prevailing wave; the 
contour line of 9-meter depth at the tip of a breakwater is same even after the completion; 
the bed slope of AC is to be 1/20, the maximum slope of deposit area of the present break- 
water. These assumptions are generally common except for the bed slope. If the entrance 
is sited near A point, shoaling will take place. To draw a figure like this the above- 
mentioned items of measurement are needed. 

In this way, the selection of the alignment of a breakwater will be soon given from 
view-point of sailing routes, tranquility within harbor and siting of wharfs. 


The way of thinking above has come from the painstaking study on relationship among 
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stal engineering, field experiments and plannings, while the author has been engaged in 
ld work. It will be fortunate if this may be a help to those engineers facing the same 
blems. As for many etual explanations omitted above, please refer to the following 


ers. 
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SHAPES OF SHORE STRUCTURES (DIKES) AND 
|'TS DAMPING EFFECT FOR WAVE FORCES 


— Especially on wave run-up and wave pressure — 


Juichi Kato 


Section of Shore Reclamation, Kyushu Agricultural Experiment Station 


This paper represents a summary of studies on rational dike design for the wave 
mping effect of various dike forms, especially on wave run-up and wave pressure. 

Up to now a number of studies have been presented on the several problems regarding 
is theme, for instance, on wave run-up, overtopping, wave forces and others. The author 
kes especially into account of the partial clapotis which varies in accordance with the 
apes of dikes and further considers on the wave run-up in company with wave pressure 
king note of the real state of locations synthetically. 

I. Laboratory Equipments and Procedures. 

Experiments were performed in a wave tank, glazed on both sides, 1 meter wide, 25 
ters long, and 1.20 meters high. The model waves generated by a ballistic pendulun, 

d the wave pressured were measures by micro-wave pressure gauge (magetic strain type) 
d were recorded by electromagnetic oscillograph with six elements at the same time con- 
nuously. 

In this study, 35 kinds of the dike models as shown in Fig. 1 that were picked out 
om about 300 sections all over the country were used and were grouped under three types. 

for the model waves, in order to generate the deep water wave and the shallow water 
ve, the angle of the pendulum swing were changed into seven kinds. That is to say, the 
lues of h,/L on Coyer’” were changed from 0.04 to 0.05. The range of shallowness, he /Lo 

from 0.00 to 0.30 and une range of steepness, Hy /Lo is 0.0054~0.072, where hy is the 
pth of water at the toe of the structure, L is the wave length, H is the wave height and 
ffix o shows the value in deep water depth. Furthermore, the wave height, lwngth and 
riod at several points in the wave tank, the height of wave run-up and spray, wave pres- 
re at the surface of dike, several dimensions of breaker, ratio of reflection and the 
her dimensions were measured. 

Moreover, the certifications of the model waves, the movements of water molecules in 
veral cases and the characteristics of breaker in case of non-structure were examined in 
tail. 

II. Wave up-rush on dikes. 


Let us show the situations of wave up-rush on dikes that were grouped into three 


pes as mentioned in the previous chapter. 
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(1) Cross Section of Dike Form, Group I (Fig.1 - 1-25). 

Cross Section of Dike Form, Group I is called the Upright Type, which is widely 
used as general shore structures in all countries. 

a) It is a well known practice to discriminate between the Clapotis Type and 
Breaking Impact Type in the incident wave on the mechanism of wave run-up on dike forn, 
Group I. The author, however, thinks that it is rational to divide further the 
Breaking Impact Type into Breaking Point Impact Type, Heighest Point Impact Type, Drop- 
ping Point Impact Type, etc. These types of incident impact mechanisr are consist- 
ently mentioned in this paper. 

b) Cross Section of Dike Form, Group I, is of various kinds because of the 
variety in the slope and the form of riprap. Therefore, the wave up-rush on dikes 
also changes correspondingly. However, it is thought that dike form may be subdivid- 
ed from the mechanism of incidence and impact of the wave into two types, i.e., the 


cross section which is composed of lines having pretty abrupt change and the cross sec- 


tion composed of smooth curves. These types give very different influence to the in- 
cident wave. Generally, on the former, the breaker is apt to occur and the spray 
flies high; on the latter, it is just the reverse. An example is shown in III, (1). 


(2) Cross Section of Dike Form, Group II (Fig.1 - 26~32). 

The dike with gentle slope, (smaller thsn about 45 degrees) on which the rate of 
reflection is small is called Group II. 

The incident situation for this type of dike is different from that for the dike 
form, Group I. That is, the advancing wave falls on the returning wave, and the ac- 
tions of plunging and knocking repeated on the sloping wall of structure uniformly. 

On the dike form, Group II, the mechanism of incident wave to the slope is constant 
even if the incident wave is clapotis or breaker type, and it may be generally present- 
ed for this group that the incident situation is regulated by the gentle slope. 

Furthermore, on the mixed form of dike, Group I and Group II, it is shown in ne 

(4). 
(3) Cross Section of Dike Form, Group III (Fig.l - 33~35). 

Dike forms belonging to Group III are made with the aim to attain the rational 
cross section df dike by the results of the above two experiments. 

The form of this Group has a large step, and the length of step and the water 
depth on the step must be determined so that dropping point of breaker may occur on the 
step. The mechanism of an incident situation for Group III should be regulated like 
in Group II, that is, double breakers occur at the end of the large step and it prog- 
resses and falls down in front of the vertical wall. Thus, the up-rush and wave 


forces on the dike wall is very small. 


After all, it is due to the difference of the rate of reflection. The situa- 
tion of up-rush and incident wave should vary with the dike form. In the following 
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(see Fig.4 - (6)) In this report, however, these relations are omitted. 

Next, the rate of wave reflection which is closely related'with the wave run-up | 
was examined. A part of the results obtained is shown in Fig.2. The rate of wave 
reflection was calculated by Healy's ceed Moreover, Fig.3, obtained by using the 
Healy's method shows the state of the wave reflection in each case. 

The wave run-up is shown according to Dr. Ishihara's ethics Furthermore, 
the relation between au/Ho and Ho/Lo with shallowness as a parameter, where ay is the 
height from still water to the maximum point of up-rush, and others were examined. 

The characteristics of the dike form were generalized. This report treats only on 
the relation between ay/Ho and hy/Lo. 
(1) Cross Section of Dike Form, Group I. 


The state of the wave run-up in each group is shown in the following three di- 


mensions: i) the limit of the overtopping of wave on a levee (sign > Auws 
4i) the limit of the fairly overtopping of sprays on a levee (sign ---- ), au,, iii) 
the limit of the highest overtopping of sprays on a levee (sign —-—- ), au7. Fur- 


thermore, B.P., i.e., the breaker point (signO) and M.P., i.e., the maximum point ( 
sign A) are expressed in every steepness, Ho/Ii9 as shown in Fig.4. 

Though the dikes in Fig.4 (1) and (2), which are very general, are very alike in 
form, the wave run-up on them presents very different values because the dike wall in 
(1) is composed of broken lines and it is interposed by curved line (2). 

To mention briefly about Fig.4, the following data on ayw/Hg and aug/Ho Or au/Ho 
are picked out: i) though spray on N - 1 (I - 2) is higher than on A (P - 1), the 
values of ayw on the former are smaller than on the latter, ii) the values of shal- 
lowness at B.P. and M.P. are very different, iii) the same as in ii) is presented 
for shallowness of the peaks in ayw- curves and aug or auy- curves, iv) the parapet 
whose cross section is composed of curves is very effective. 

Moreover, the curves (sign--- -- ) shown within Fig.4 - A that were calculated by 
, Stoker's impact wave theory, are not a little different from the curves of the wave run- 
up on the dike A. 

As for the dike under the construction, its stability for the wave pressure was 
investigated using three kinds of the reduced scale, 1/10~1/30. ~ But it is omitted 
here. 

(2) Cross Section of Dike Form, Group II. 

The curves in Fig.4 (3~4) show the wave run-up on various kinds of 1 : 2.5 

slope which belongs to gentle slope type, where Azs- s is a dike with a sigle slope 1 : 
2.5 and Ags- a, Azs—- b and Aps— c are dikes with vertical wall of same level height on 
the slope 1 : 2.5 of Aps- s dike at the different positions respectively. 


Wand $) 
As for Ae5- s the same results as by Dr. Sat6 were obtained. The value of 
a@ shows the maximum point of the wave run-up. How does the dike form with the ver- 
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Fig3 RATES OF WAVE REFLECTION 
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tical wall on the slope influence the wave run-up ? Now let us compare this result 
with those in Aas - a, Ags - b and Azg5- c, etc. ; | 

The data obtained from these curves are the following: i) the maximum peak of 
curves occurs on the curve of Hg/Lo = 0.04, ii) as the value of shallowness becomes 
greater, the influence of the upright wall appears, iii) as the upright wall moves 
backward such as Aps - b, Aas - c, aug/Ho curves tend to approach the run-up in Aas - 8. 

So it may be assumed that the curve in Ag5 - s corresponds with aus/Ho curve in 
other cross sections. Thereupon the auw/Ho curves in the Apc - s cross section were 
decided according to the appearances of the plunging waves on the slopes. Consequent-— 
ly, it was judged that these curves are approximarely equal to the ayw/Ho curves in Azz 
-c. In each graph the curves in Az; - s are shown with smaller broken lines for 
reference. 

Here, special mention is given to the fact that as the rate of reflection be- 
comes larger, B.P. occurs after M.P. (the same is true in Group III). 

From these results, many useful data in practice can be obtained. For in- 
stance, it is desirable that the form of cross section with broken lines, such as a 
berm, is avoided. 

(3) Cross Section of Dike Form, Group III. 

As for this type, even the largest value of ay7y/Ho is only about three. The 
dike in this Group has a pretty similar character to that of Group II. Furthermore, 
as if there is little difference between the cross section R- 1- a and R-1 = b. 


The arrow marks in R - 1 graph show an example of the height of spray, 7% , when 
breakers fall on the dike wall in case of large h and the breaker impacted on the up- 
right wall. As in this case the appearance of the run-up is like that of Group I, 
the property of Group III vanishes. The length of breaker zone, 1, can be obtained 
from the same graph. And the arrow marks in R —- 2 show the height of spray,7’, when 
the breaker impacted on the slope of the second step. The selection of the stair 
having one step or two steps must be decided from the construction for the foundation 
and the waves. 

(4) General Cross Section of the Dike Form ( On Use of Single Slope Method ). 

In regard to the representative dikes as mentioned above, the wave run-up was 
analized from the data of the experiment. But, how can this problem be solved con- 
sidering the general dike containing the dikes in Group I, II, and III ? It is con- 
venient if the height of the wave run-up can be estimated by other methods. 

On this problen, sativieienepedes an interesting result by the Single Slope 
Method - " when incident wave becomes a breaker, the height of wave run-up on general 
composed structures is equal to the height of wave run-up on a single slope which is 
drawn from the breaking point to the maximum point of wave run-up. " 


So, the author examined the results of the experiment by this method and the re- 
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ts were as expected. Here, the author expanded the method and arrived at the fol-— 
ing conclusion: * 

"The height of wave run-up on dike can be replaced by the height of the wave 
iup on a single slope which is drawn from the breaking point in Ber of the break- 
; incident wave or from the point of L/2 length from the dike in case of clapotis 
ident wave, to the maximum point of wave run-up. " 

Fig.5 shows the up-rush on a single slope. Then, according to Fig.5, using 
» Single Slope Method, the up-rush on optional dikes such as Group I, II, and III, 


;e, Can be estimated. Fig.6 obtained by the same method for dike form as shown 


yve illustrates an example showing its accuracy. 


IV. Intensity of the Wave Forces. 

To study on the effect of wave damping for several dike forms, the intensity of 
re forces should also be considered taking into account the wave run-up simultane- 
sly. Many people have studied on the wave force, nevertheless, up to now, there 
> but a few practical methods available. On the wave pressures, the theory of 
‘tial clapotis is not established, so the important relations remain obscure as the 
_lowing: i) How is the problem of wave reflection related to the wave pressure ? 
.) How is the wave pressure of its medium type in incidence between clapotis and 
aker which takes place frequently in the actual wave be calculated ? iii) How is 
> relation between the various mechanisms and the wave pressure treated ? In this 
ort the rational formula based upon the numerous results of the experiment is in- 
ed. 

(1) Results Obtained. 

The following relations were considered from the results of wave experiments: 

' The vertical msximum wave pressure distribution. ii) Relation between Bro/Wo 
|Ho.  ii1) Pmox/Wo and Ho/Lo. iv) Bnav/Wo and C°/2¢. V) Proy/Wo and T (where Prox 
maximum wave pressure, C is the velocity of the wave, T is the period of the wave). 
re, only remarkable data are enumerated. They are: 

On the breaker: i) In this experiment, there was no occurrence of the phe- 
enon reported by Manik? that the crest of the wave engulfs tne air and dashes a- 
nst the dike wall compressing the air horizintally. But it is considered that 

phenomenon, for the crest of the wave to engulf, may occur occasionally only when 

ure progressive wave with small value of steepness dashes agsinst Group I dike 

se riprep is composed of broken line in its cross section. So it is thought that 
occurrence of such phenomenon as Minikin reported is possible actually by the in- 

ference of the weves themselves. ii) The actual wove pressure is accompanied by 

oscillation owing to the complex reflected wave. iii) The maximum pressure dis- 
bution dixgrams in Group II and III have a definite trend clearer than in Group I. 

) The relations between Pmaxand Ho, Ho/Lo and C*/26 have a definite tendency on the 
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nd of incident mechanism respectively. These results up to now are not very 

ir. v) Both the positive and the negative force in the impact on the dike step 
‘Group II and III interact frequently at the same time. Therefore, the attention 
have to be paid to the construction of dike. 

| Next, on the clapotis: vi) With regard to the vertical maximum pressure dis- 
pution for Group I, the maximum value of wave pressure is larger than the theoreti- 


value computed by Sainflou's formula. 


After all, it can be concluded that the typical continuous or non-continuous or 
1lsive periodical fluctuation as reported by Varyaee or by seineion? does not occur 
weve pressure diagran. 

(2) The Formula for the Maximum Pressure of Breaker on Dikes. 

As the mechanism of breakers is very complicated, it is difficult to solve prob- 
s on the breakers directly from the standpoint of hydrodynamics. However, accord-— 
to numerous pressure data obtained by this experiment, all situations of breaker 
act may be able to consider as the phenomenon that a jet dashes, in spite of the 
lous incident situations (i.e., in spite of the various dike forms). So, con- 
sring this phenomenon as a jet flow in the same way as Dr. Nagai's noribae the 
Pficient of the maximum pressure formula is determined experimentally in each im- 
>; Situation of every cross section of dike: 

That is 
os aa pce Where p is density of water. 

This formula shows the pressure of a jet flow which dashes against a vertical 
DB However, as the motion of the water particle which varies according to the 
> form is very complex, Pmaxis shown by the equation containing the coefficient, wu 


srmined in accordance to each dike form. That is 


Ze P max = P Co IVS bm 


where Co is constant. 


Now, the first approximate value of the theory of shallow water is shown as 


lows 


3. Ra= tankhkh 


2 


where Q is the half wave height and Ris TZ 


A. lL = 2NQA/tanhkke 
If the wave length at the water depth Ai [= f(Hv)] where the incident wave is 
isformed (mainly, it is the water depth at the toe of the dike) is shown with Ly}, 


coefficient of dike form is given by the formula 


Y ee 
De = See 
aE AES aa 
Then, if Coit, is replaced by Ll, Praxis given by the formula 


6. Pmax / Wo i MT ApHp / Lb where ic is VGhb ana 2a LS equel to Hp (hy). 
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This formula is the general form of the computational equation for intensity of 
the maximum pressure of breakers. In Fig.7 the incident situations are classified 
in each group of dike according to the experimental data, and the relation between Inox 
and h,H,/L, is shown. The values of u“ can be decided from Fig.7. Each line in 
Fig.7? is shown according to the classification of the incident situations correspond- 
ing with the range of steepness to each shallowness which is obvious in Fig.4. 

Furthermore, according to this experiment, the distribution of wave pressure in 


each case can be decided by using equation 6. as shown in Fig.8. 


V. Rational Cross Section of Dike. 

From the above-mentioned, we know that the wave action on dikes varies with the 
remarkable change of condition in the locations. For instance, first of all, the 
depth of water at the toe of the structure, wave shape, interference of wave itself, — 
and moreover, wind force, etc., and that at the some time, these elements are affected 
by the various forms of dike cross section. So, we have to design dikes taking into 
consideration that the form of dike must be selected so that it can be expressed in 
the most advantageous terms. 

That is, the wave energy must be dissipated not by shore structures directly, 
but by the transformation of wave according to the shape of dike, i.e., by making the 
breaker waste its energy on account of the turbulence. So, it is necessary to make 
the dike in such a way that its cross section promotes the formation of breaker and 
its turbulent energy exhausted just in front of the dike wall. Moreover, the shape 


of the dike is influenced largely by the condition of the location, for instance, geo- 


graphical features, foundation of site, etc. So, we have to decide the rational 
cross section of dike from all angles. Thus, we are made to realize that not only 
one rational shape of dike can be brought into existence in the future. Cross sec-— 


tions of dikes from Group I to III have merits 1s well as demerits, and we must choose 


the most profitable shape of dike according to circumstances. 

From this, we can gain three types of the rational cross sections of dike as 
shown in Fig.9. 

Let us meke note on Fig.9: i) The surface wall of the dike in (A), (B) and 
(C) type is constructed smoothly by curves so that the slope composed of riprap at the 
dike's toe may not change abruptly, and the dike itself may do the duty of a parapet 
in order to control spray as much as possible. ii) The three dike forms (A), (B), 
and (C) are successively lower in height and longer in width in case of constant depth 
of water at the toe of the dike. Fig.9 is shown as the example with water depth of 


3.6 meters at the toe of the structure, the wave height of 2 meters, the wave length 


of 50 meters and steepness of 0.04. iii) The characteristic of type (A) is that the 


cost of embankment is the cheapest “mong the three, and the experience of construction 


is abundant now-a-days, although in ceneral, its wave run-up snd wave pressure are the 
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2st. So the banking of type (A) is unsuitable in the site where the foundation 
oft and the change of water level is very large and where typhoons often occur. 
the other hand, type (B) and (C) are very suitable for weve forces, but the cost 
nbankment is not at all cheap. So these two types are unfit in ordinary founda- 
and where the change of water level is small. iv) The problem on the general 
truction, the form of rear part of dike and the others (including the problem of 
2ability and drainage, etc.) is omitted in this paper. 

The author would like to reiterate that this is only an example for the princi- 
9f dike form and the dimensions of these types are changeable adequately in accord- 


with the circumstances. 


In this paper, the author shows the characters of incident wave to shore struc- 
s and of wave run-up and wave pressure in every case. Here we can derive a prin- 
> of practical dike design. As mentioned in the preface, this is only a summary 
ne study, so the author would like that Reference 11) in which this theme is dis- 
ed in detail be read. Moreover, in this study, each problem is considered only 
ne situation of no wind at two dimensions, so the studies in which the case of 
gue incidence and problem of diffraction and reflection of wave and wind forces, 


will be taken into account are expected in the future. 
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ynopsis: This report concerns with the estimation of the shore line deformation of sand 
eaches by the variation of longshore current velocity along the shore line, when the 
ffshore breakwater is constructed. The development and final stable forms of tombolos 
esulting from the construction of offshore breakwater are experimentally analyzed in 
erms of the influences of length of breakwater, water depth at which the breakwater is 
onstructed, anc steepness of deep water waves. 

Furthermore referring to these basic experimental results, the deformation of shore 
ines at Miyazu Coast, kyoto Prefecture, is estimated, when the offshore breakwater is 


lanned. 


- Introducticn 

Among many tyves of shore deforriations in beach processes caused by sand drifts, the 
each deformation in short duraticn results mainly from the movement of sand drifts in the 
frecticon of beach profile, whereas that in long duraticn such as a year, several years or 
ore, is com only influenced by the littoral sand drifts with some exception in several 
ases of specific coastatal features. When coastal structures like jetties, breakwaters, 
ffshore breakwaters and sea dikes are planned, most significant is that the future shore 
satures after completion of coastal structures must be estimated es possible as accuretely 
> insure the control of the beach deformation against the wave xticn. i’roblems, therefcre, 
sncerr.ed with the effects of coastal structures on the littoral sand drifts beccme important 
> secure the hydraulic performances in beach protecticn, though each cf ccastal structures 
s distinctively different. 

Many Meee ve been made for the effects of shore structures on the defcrmation 
> shore line. These studies, however, treat mostly with the performances of jetties, beach 
soins, or levees on the beach process, whereas the influences of offshore breakwaters and 


36 dikes on the beach deformation have scarcely been studied. Offshore breakwaters are 
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frequently used as a mean for prctecting anchorages and beaches, and their effects on the 
littoral sanc drifts are of quite irpcertance. 

For clarification of the effect of offshore breakwater on send drifts, the following 
subjects must be studied: (1) duration and amount of littoral sand drifts, and (2) zone 
of movement of drifts. The relation between the amount of littoral sanc drifts and the 
velocity of longshore current, which is the most imoortant factor to the moverent of sand 
drifts, has not been established, anc therefore nast studies made by the authors were only 
to estimate qualitatively the deformation of coastal features from the velocity distribu- 
tion of longshore currents. 


In this renort, first is ascertained from the exneriments conducted in the model basin 


that the method of estimation of the deformation nrocess in the ccestal features by the velo 
city distribution of longshore current is available for engineering needs and an anelysis of 
the effects of offshore breakwater length, water cepth at the bre-alwater and steepness of 
deep water wave on the formation of tombolos and its stable shaves is presented. Further- 
more, with results of exverimental data and the littoral distribution of the amounts of lit- 
toral sand drifts commuted fror wave energy, the deformation process at Miyazu Cost and its 


final feature are estimated, when the offshore breakwaters are constructed. 


2. Transformation of wave heights by refraction and diffraction 

When the waves come obliquely from offshore with a definite angle, they are sradually 
refracted by the topographical features. If there exists an offshore breakweter near the 
coast, the refracted waves are also influenced by the diffraction due to the structure and 
thus the wave height, the incidental anrle and the bresking depth in sheltered area are ap- 
preciably changed. The sare diffracted behaviours in wave characters are evidently seen, 
when the waves come normally to the structures. The estimation of changes in wave chsrac- 
teristics by refraction and diffraction is of basic importance to calculate the wave energy | 
and the velocity of longshore current in the direction of shore line, therefore the change 
in characteristics of breaker refracted and diffracted by offshore breakyaters must be re- 
vealed for clear estimation of beach deformations in sheltered area. Letting the equivalent 
wave height to deep water wave (not refracted and diffracted) denote by H'y the following 


expressions are obtained. 


H-+(2) Hp =f (te) (1) 
H, L.o/, H, i 

where H is the wave height at the denth h, L the length of deen water wave, Hp the bresker 
height and hp the breaking depth. 


The relationship between H and Hp, is expressed as follows, when the incom/ng waves are 


¢ 
| 
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racted and diffracted, 
tsi = bh Hp = ; 
asain Mb, te) 


re K is a coefficient of refraction at water depth h, and K'p a coefficient of diffrac- 


n at breaking point. Tne ratio of Hp to Hp is then, 


Lh an Vay Ga 0 mo VQ h Dp ( 

He = KK) = KK CE) / FES) J, 3) 
inserting Eq. (1) into Eq. (3), tne eauivalent wave height to deep water is expressed 

/ 


KK, {4(2) /4 8) = 


"in Eqs. (3) and (4) are evidently functions of water depth, Hp/Hp can not 


= 


ce K K'p and f 
obtained only when H,/Lo is given. The relations of H)/L, and hp/Lg to K K'y for partic- 
r waves of definite period must be evaluated in the region where the wves are diffracted, 


can be calculated by the following procedures. 


(1) Assume the wave period T or deep water wave length Lo, and breaking water depth hp. 

(2) Obtain H'o/Lo, Ho/H'o f'(hp/Lo) for particular waves and particular breaking depths, 
and H/H' )=f(h/Lo) for the water depth at both ends of offshore breakwater. 

(3) Compute the coefficients of refraction and diffraction for hy assumed in (1), by 
the following method. 

(4) Using the values obtained in (2) and (3), compute the steepness of deep water wave 


by the following equation of (5) for given values of hp, 


ihe Walp Se | (be) / #2) } (5) 
ee SlsoekK Ky Lio! bso 
The coefficients of rofraction and diffraction in (3) are obtained with the use of 
raction and ditfraction diagrams. If the beach is constant grade and thus all contours 
parallel, the refraction coefficient can be obtained by the law of Snell, without using 
refraction diagram. Denoting the angle between the crest line of waves and the contour 
e at the water depth h by o , and the angle between the crest line of deep water wave 


the contour line by; oo? the following relations are obtained; 


Sih eae 27h (6) 
Sint camen Cs lam h one 

/ 

Daal ay WM 
K alee) i ian as 


re § is the distance between two orthogonals of wave crest line at waterdepth h, and 
fix o indicates the one for deep water. With the use of these equations, K is cal- 
ated for definite values cf d,and h/by at all points. frig. 1 indicates the relation- 


3 6) 
p of these three variables, &); n’bo and kK. or tue diffraction, J. W. Johnson has 
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calculated the changes of wave heights and phases from three to 360 wave lengths and J.W. 
Dunham also from one to 6 wave lengths, when the waves come perpendicularly to an infinite 
peninsula. 

Offshore breakwaters are commonly con- 
structed within the distance of several wave 


length from coast, so that the detailed compu- 


tation for wave diffraction within the distance 


of two wave lengths was made by the use of 


8) 
Sommerfeld's exact solution. 


Fig. 2 is the diffraction diagram obtained 


by the calculated results, and also the experi- 


mental ones in wave height in a form of the ast 


ratio between the wave height in the sheltered 


(oy yea a en 
area and that’ at the offshore side of break- 56 8 Q0/ 2 Sb Be: (0) 
water are indicated by numerical values near Mh 
small circles. As seen in this figure, a good 
agreement between computed and experimental Fig, 1 Relationship between «, , 


h/L, and K 
values for wave height is obtained. 
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Fic, 2 Refraction diagran in the case of y/L < 2,0 
3. Beaches deformation in sheltered area and formation and development of tombolos 
: 
The veiocity of longshore current can be obtained by the height of breaker and the 


incident angle. In the past, beaches in sheltered area were presumed to be eroded if the 
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plocity of longshore current increases in the direction of running current and deposited 
| 


- the velocity decreases. 

In this chapter, the validity of the preceding method of estimation for beach deforma- 
ton by the velccity distribution of longshore currents will be concerned, and the formaticn 
r tombolos and its develcpment as results of beach deformations are also investigated. 

1) Velocity distribution cf longshore current and deformation of beaches 


in sheltered area. 


The model besin of fun shape shown in Fig. 3 was used for vericus runs of experiments 


© that incident angles of 100 
aves can be changed. The 

SS 
elaticnship between the ee 
elocity distribution of & 

™ 60 
ongshore current and the & 

Ss 
eformation of beaches was ue 40 

Ww) 
xperimentally stvdied, when e 
he offshore breakwater was 4s 20 

pS 
laced in the straight shore 


f constant bottom grede. 


O 
Ql 0203 060810 20 
Sieve Opening mm 
Fig, 3 Flane of the model Fig, 4 Sieve analysis curve 
ype were used for velocity basin of fun shape of used sands 


urrent meters of propellar 


easurements of longshore current anc the records 
ere made by the electrical devices. 
= 40) = 00497 


The mean diameter of sands used is 0.42 mn, ony s 
c (nl. 8-0 @ 


he sorting coefficient Jaq5 7425 46 1.26, the 8 400 
pecific weight is 2.66, and the grains of sands Scot 
re approximately uniform having the grain size Q 

f sieve anelysis curve as seen in Fig. 4. s z 
xperimental waves generated with angles of 30. s 
nd 4O from shore line. The offshore breakwaters ® a 


ere placed at angles (©) of O°, 30° and 40° from 
hore line for all waves. 


Fig. 5 indicates the relation between tne 


Veloci ty of 


Longshore Curren 
S 


' 
Ww 


mount of sand transported per unit width in va- 


j heltered area and the velocity ’ : ‘ 
Bae there trene 5 Velocity distribution of 


longshore current and ex- 
perimental results of beach 
It is process in the case of = 40° 


istribution of lonyshore current along a line 20 
mn apart from the original shore line. 


vident that the wave height, decreéses 
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gradually in the inner part of sheltered area with results of decreases in longshore current 
velocity as estimated by the description in (2). In Fig. 5 wnich indicate the case of X= LO” 
positive and negative signs in the velocity of longshore current indicates the components 
along the shore line in the progressive direction of waves and in the reverse direction, 
respectively, and O denotes the angle between the offshore breakwater and the shore line. It 
is seen in this figure that the positive slope in curves of the velocity of longshore current 
indicates the possibility of erosion, and converse'y the negative slope indicates that of 
deposition in both cases of 9=0° and 40? The ercded area and deposited area, however, are 
largely influenced by the angle between breakwater and shore line. For 9=0O, the erosion oc- 


curs at the frontal side of breakwater facing the progression of waves and the maximu:. 


deposition is seen near the other side. On the contrary, 
when the breakwater is placed at 6@=40 and parallel to the 
crest line of incident waves, the maximum deposition occurs 
at the former side, whereas the erosion is seen at the lat- 
ter side, and it indicates the opposite tendency to the 
case of @=0. 


It will be, therefore, evident that the offshore break- 


waters constructed at a definite angle to the shore line and 


resulted in the beach deformation, give much influence to 


the final tongue-like shape of tombolo which is a result of 


Fig, 6 Representation of 
beach progress due to erosion and deposition. It will be symbols used in 
the study 


concluded that the future beach deformation can be forecasted 


for engineering purpose, if the velocity distribution of longshore current is estimated. 


(2) On the formation and development of tombolos. 

The formation and development of tombolos have been briefly described in (1), by means 
of the diffraction phenomena of waves. In this section, influences of the steepness of deep 
water wave, the length of offshore breakwater, the water depth at the location of breakwater 
to the formation and development of tombolos are concerned, and furthermore, the final estab- 
lishment of tomboloe elso is treated by the experimental results for A= O-and @= eee 

i) Influence of length of offshcre breakwaters: 

Before the analysis is concerned, symbols in the study are explained in Fig. 6, h, is 
the water depth at the location of offshore breakwater, and Hy the wave height at the front 
side of offshore breakwater. The influence of the length of offshore breakwater was in- 
vestigated in the four cases of B/Lp=3.00, 1.00, 0.50 and ©.25 under the definite condition 
of Hp/Lp=0.05, h/l=0.C5 and T=0.807 sec. The experimental runs revealed that the change 


of amount Vi, of deposited sand in the sheltered area behind the breakwater with time. 
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ii) Influence of location of offshore breakwater: 
| As seen in diffraction diagram of Chapter 2, the equivalent deep water wave steepness 
aries gradually with distance between the offshore breakwater and the beach line. The 
ave height H after diffraction is given by 

H = K'ces2x (4 -/ ) (8) 

H, is 5) : 
here K' denotes the coefficient of diffraction at water depth Cys phase coefficient, 
the time, and T the period of wave. 

The rate of sand drift by longshore currents parallel to the shore line at x/L, from 
he offshore breakwater is empirically ex- 
ressed by the following equation derived 
rom the field observation at California 
Bein” end Miyazu coasts as well as ex- 


erimental results. 


“™m ; 
ty = RE; 59) 
here Ej is the energy of waves in the 


irection of shore line and # and m are 


onstant. Since #4; and cos Y are, 


EF oc H*L,coo¥, 
COI YU oc Cop) 


(10) 


here denotes the angle between the 
nore line and the progressive direction 
Fs ODS, . 
f diffracted wave. ‘ir, 7 Relation between k "(-M4+) 
y/ly and gli 


when the incoming waves have the same 


characteristics, Eq. (9) is rewritten «cs follows, 


m 
( 
Gy oc Pa ON) e (11) 
ted ay 
p) 

Curves of K'? Cay) to y/Lo are indicated as a parameteric expression of x/Lo in 
ig. 7. Curves are steep for small values of x/Lo and tney become gradually flat with the 
nerease of x/Lo. Since qy is related with the x? Cafd/y)turcugh Bq. (11), the heavy rates 
fF sand movements are indicated by steep curves in Fig. 7, and therefore fcr small values 
fF x/Lo, the transportations of littoral sand drifts, which results in erosion and deposi- 
ion, are appreciable. 

No large variation in transportation of sand drifts along the beach line is prceduced 


y the longshore current, despite of the movement of sand particles, when the offshore break- 
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water is placed far from the shore line, and the sand drifts are little deposited in the 
sheltered area. On the contrary, the sand drifts transported decrease rapidly with the 
increase of y/Lg in the sheltered area when the offshore breakwater is placed near the 


shore line, with results of appreciable deposition. 


iii) Influence of steepness of deep water 100 
wave: 

The wave of Hj/Lo = 0.05 in this experimental 80 
run is transformed into the wave of equivalent deep 
water wave steepness H')/lo< 0.03, at which the sands 4 
are extremely transported. The wave of H,/Lo= 0.02 %o 
becomes the equivalent deep water wave of H',/lo < a 
0.012 gives less movement of sands compared with the 
former case. Therefore it is estimated the rate of = 
deposition is larger in the case of H,/L,= 0.05 than i 


in that of 0.02, despite of the value of F,/Lo. 


This facts is revealed by Fig. 8, which shows the 


rate of deposition at various cross sections in Fig, 8 Variation of rate of de- 
xi positicn by chanee in 
sheltered area with the changes of Ho/Lo and ho/Lo. Hof, and ho/L, 


iv) Final shape: 

the formation and development of tombolos 
have been discussed in i) ~ iii), and in this 
section, the final shape after establishment 


4 ec 
of tombolos will be treated. Fig. 9 shows the pease 


offshore tle = 005,P"-a05 
Breakwater 


final shapes of the shore line for various 
values B/Lg and it is seen that the resulting 
beach line is indevendent of the length of off- 


shore breakwater. When the depth at the loca- 


tion of offshore breakwater is taken as vari- : ‘ 
Bip, 2 Final shapes of shore line 


able, the final features of shore line are of for variable 3/ig 


similar shapes, independing on the shallowness, h./L,+ The shape of shore line coincides 

to the crest line of incoming waves after diffraction. Therefore, it is considered the 
final shape of tombolos will be estimated, when the incident waves come normal to the shore 
line and no sand drifts along the snore line are produced. Since the shane of tombolos be- 
comes stable when the shore line coincide with the crest line of diffracted waves, there 

will be some relationship between ‘he dentn of water ho at the location of offshore break- 
water and the distance B' from the intersecting point of tombolo and offshore breakwater to 


the end of structure. The experimental relationship is indicated in Fig. 10. Thus, with 
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e use of the diffraction diagram, the future a! 


x10 

atures of shore line will be estimated for Mw 

ven water waves L,. 

| The influence of H,/Ly to the final shape B 

. tombolos will be next considered. The Te 

Baas of diffraction diagram are not in- : 

Luenced by the values of H,/L,, as explained 4 

y the theoretical solution on diffraction by 

ommerfeld. Although the values of H,/L, in- | 

luence to the rate of growth of tombolos, its 

inal features will be considered not to be O < 
a) 2 yi 6 8 jor? 

afluenced by the wave steepness of deep water he/p . 

AVES. 

Fig. 11 shows changes in the rate of ed ae eee Bile 


=position at several cross sections in the 
neltered area with the change in H,/L, and 


ne evidence of foregoing description, and ; 
>» significant influence of H,/L, to the 


© 
o + -0050,18h" 


inal shape of tombolos will be seen. 


The beach deformation in the sheltered 
rea by the oction of diffracted wave for 
.= 6 and @=0 has been discussed in above 


sections. Problems on the final shapes for 


* C,@6 * O remain unsolved. Further 
asic studies should be made in order to 


larify the process of beach deformations in 


ris case. 
. Bstimation of beach process at Miyazu ‘ig. 12 Chan-es in rate of dénosi- 
tion et each cross secticn 
Coast net oe SA 1S te Al . 
Ur 1e) 


In this chapter, using the foregoing 
asic experimental results, the process of beach deformation by the construction of off- 
nore breakwater at Eziri of Miyazu Coast, Myoto Prefecture, is revealed. «as the first 
teo, the authors estimate the beac’ deformation from the velocity distribution of long- 
nore current along shore line and check the estimation by the distribution of amount of 


and drifts, along shore line, using the relation between the wave energy and the amount 


> sand drifts, obtained by the field observations. 
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Furthermore, the final shapes of shore line in sheltéred area resulting from the con- 
struction of offshore breakwater will be treated, using the descriptions in Fig. 10 and the 
diagram of diffraction. 

(1) Estimation of beach deformation by the velocity distribution of longshore current. 

Miyazu Coast situates in the west coast of Miyazu Bay, the mouth of which opens north- 
east toward the Japan Sea, and Amanohashidate which is a typical sand spit is in this coast. 
The bottom grade is nearly constant. The sand drifts in this coast are transported by waves 
of incident angle 85° to the shore line ané moved from north to south. 


The period of incident wave at Miyazu Coast 


i timated sec by the analysis cf wind 
Ware drecton XX Aah ies 37 y past ; 


Break water velocities and directions. 

The offsnore breakwater for anchorage shown 
in tig. 12 has been planned and constructed 
partly. Before the velocity of longshore current 
in the sheltered area by the offshore breakwater 


is estimated, the variation of breaker height 


and the breaking angle to the shore line result- 
ing from tne refraction and diffraction must be 
“ig, 12 Location of offshore break- ? 

weter at Miyazu Coast definitely evaluated. Fig. 13 is the relation- 
ship of H,/Lo, hy/Ly to K*K', for Miyazu Coast 
described in chapter 2 and it is seen in the 


figure that the breaker height and breaking depth 


for a particular deep water wave cf a definite 


steepness are variables of locations in the shel- — 


tered area. In the figures, the values of y re- 


present the distances along the shore line as in 


icisleaey aes 


The coefficient of friction ky must be known 


for computation of longsnore velocity. Inman and 


1D 
Quinn” proposed the empirical equation of (12), 


using a large number of data obtained by field 


re pekaraea e. aE ERNE observations and experiments, 
Ppl ante ihe wang. Ullllee 
[Om Zan TEES Oem eek CRe!/O ‘ a ap, F 
hey kg = 0.00401 ‘ Cm = see 2 (12) 
= 0.024 vy # Cit. see), 
Fig. 13 Relation H./L,, ho/L, an: 
KK! at Miyazu Coast in The calculated velocity by uc. (12), however, is 


the case of T 3 sec. 
different from the observed results by the authors 
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that the authors obtained the value of kg by the 
Heerved data through the equation of Putnam, Munk 


12) 
mad Traylor on the velocity of longshore current. 


ig. 14 shows the relationship between kf and ob- 
erved velocity of longshore current and the follow 


ng empirical equation for kr will be proposed. 


kg = 0.0249 Vv (m= sec ) 


(fth=ssec))) (13) ome 


NIG NIG 


= 0.1476 V- 


lien 


. 2. 
ulated for particular waves of definite characteris- i Ue OA GeGsi0" 20" 40 G0G00 
Observed Velocity HKec. 


sing Eq. (13), the velocity of longshore current 


t various locations in the sheltered area are cal- 


ics in period and steepness, and the results are 


hown in Fig. 15. The regions of erosion and deposi- Fig, 14 Relation between ccef- 
ficient of friction end 
ion estimated by the decrease and increase of long- velocity of longshore 


hore current velocities are also indicated in the same figure. 

As seen in this figure, when an offshore breakwater is constructed at sandy beach, the 
and drift is deposited in the inner part of breakwater, whereas the erosion is occurred in 
he out side of breakwater, independing on the characteristics of incoming waves. 

2) Estimation of beach deformation by the distribution of amount of littoral 

sand drifts. 

Fig. 16 indicates the relation between une transmitted energy of breaker per unit length 
iong the beach line expressed by Eq. (14) and the total amount of littoral sand drifts Qy 


n terms of a parametric expression of H')/lo. 


E. = -wHy(Ls/T)sin2a,, (143 


1 
cm/s 
[00 
where Lp is the wave length of breaker, 
80 Te75 t2=002 
Sie. Spore w the specific weight of water. Fig. 17 
ie \ T=6, (25002 shows the distribution of amount of lit- 
40 toral sand drifts for the wave of T= 5 
20 sec and Ho/Lo = 0.03, and that of T = 6 
O sec and H)/L, = v.02 with use of K*K'y 
O GO 120 /80 240 300 


obtained in (1), By in Eq. (14) and Fig. 


Distance along shore Line (M) 


Deposition | P Erosion 28 16. Evidently seen in this figure the 


prediction of beach features by the 


Mig, 15 Distribution of longshoere current Gistributicon of amount of littoral sand 
velocity in sheltering erea 


drifts is reduced to the same result 
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estimated by the velocity distribution 


Q ar 
of longshore current. A 2 
Z O Ex hve od =10 
+ w «20° 
(3) Final shape y 3 30 
e@ . a 240" 
As shown in Fig. 12, the breakwater Oipetteb tinh: 


2) Suffix Numbers alate Spas 


is made up with two parts; the one @ 


runs towards the shore and the another 


@) runs parallel with the shore line. 


As regards the part @) of the offshore 


breakwater, resulting function of the 


tombolo to the littoral sand drifts be- a 
| ; 
comes the same one as of groin. Then Bee cas IL | a | 
Se f , 
we have only to consider the diffracted ae Z 7G re ein ae Pi me COR 4 6 8/0 
Ei: Sm 


wave by the offshore breakwater @) 

which can also form a tombolo. The final 

shape of this shore will be able to be 40 

estimated, since the final shape of 

tombolo is not influenced by Ho/lo and 

B/Lo as described in the previous chapter mM) 

and is determined only by h,/L,- é 
In the plan on construction works of ZO 

breakwaters at Miyazu Coast, the depth of 

breakwater placement is designed to be 10 

3.5m. The final shape of tombolos there- 


fore is determined only by the wave length 


0 
L.. Table 1 indicates the value of B' for O 60 120 /80 240 300 
Distance along shore Line (m) 
Deposition el Erosion | Dep. 
Fig, 17 Distribution of emount of littoral 
send drifts in sheltered area 
! 
Fie, 18 Final shape of shore line in Table 1 Values of B 


the case of T = 6 sec 
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ives of 5~7 sec. obtained fiim Fig. 10. The value of B' becomes large with decrease 


e length as seen in Table 1. and thus the tombolos hardly when wave length becomes 


rig. 18 shows the final shave of tombolc for the wave of T=6 sec. which is most commont 


yazu Coast. When the offshore breakwater is designed as anchorage, if the depth at the 
nent of offshore breakwater is constant, the anchored area is determined by the wave 

1, so that longer breakwaters than those determined by the previous discussion are un- 

sary for coastal design. 

is seen in Fig. 18, a good agreement between the final shape predicted by the empirical 
ionship described in Fig. 11 and the estimated result by the distribution of longshore 

1t velocity and the distribution of amount of littoral sand drifts is obtained. 


ynclusion 

In this study, the authors proposed the method of estimation for beach deformations by 
-fshore breakwater by means of the velocity distribution of longshore current, and it 
srified from the experimental study and the field observations of sand drifts at the 

, The effects of the breaker length, the water depth at the location of breaker and 
beepness of deep water waves on the formation and development of tombolcs were experi- 
ily investigated, and the final stable shapes of tombolos were also treated. With the 
> the results obtained, the process in beach deformation at Miyazu Coast was estimated, 
the offshore breakwater would be constructed. 

further investigation to obtain the influences of incident wave angle and the location 
*shore breakwater to the beach deformation by the basic experiments will be extensively 
ipate. 

The authors desire to express their hearty thanks to Prof. Ishihara and Assist. Prof. 


ci for very valuable advices and suggestions during the work. 
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_ ARRANGEMENT OF GROINS ON A SANDY BEACH 


Shoshichiro Nagai” 


It is an important and difficult problem to arrange groins effectively for protection 
st erosion ty wave action on a sandy coast. This paper will present the relation of 

roins' length, space, and orientation with respect to the shoreline, the direction of 

propagation, and the breaking point of the breakers. The relationship between wave 

ness and sand transport, and some results of experiments concerning special types of 
s are also presented. 

Groins, constructed perpendicular to, or at an inclination with the shoreline, have 

used as one of the most effective measures for coastal protection on a shore which 

n abundant littoral drift coast. Notwithstanding a few laboratory experiments and 
investigations with the groins, it has until today been impossible to deal rational- 
th the length, direction with respect to wave climate, space and arrangement of the 

$s on sea shores. Therefore only practical experience has played an important role in 

problen. 

The writer has been studying since 1954 about the relationship between the length of 
s, the space between them, the orientation with respect to wave propagation, the type 

haracteristics of waves, and the longshore currents due to wave breaking in the Field 

ulic Laboratory of Osaka City University. Only the results of the experiments on 


s will be described in this paper. 


MOTION OF SAND PARTICLES BETWEEN GROINS 


Shoshichiro Nagai and Hirokazu Kubo** 


Experiments on the motion of sand particles between groins were performed in a fixed 
to compare with the results of groins in a movable bed. The comparison proved that 
results were in comparatively good agreement. 


Arrangements of groins on a sandy beach were studied for two years, 1954 and 1955, 


sf. of River and Harbor Engr., Faculty of Engr., Osaka City University. 
st. of River and Harbor Engr., Faculty of Engr., Osaka City University. 
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in a basin with a movable bed, and the results .ere reported in the Proceedings of the 
American Society of Civil Engineers, Separate Paper 1063, WW 4, September 1956. The 
experiments dealing with the motion of sand particles between groins were made in a basin 
with a fixed bed during the period April 1957 to March 1953. This paper presents a com- 
parison of the results of the previous movable bed experiments with the results obtained 


in the fixed bed experiments. 


THE SHOCK PRESSURES EXERTED BY BREAKING WAVES 
ON BREAKWATERS 


Shoshichiro Nagai 


Since the prediction of shock-pressures exerted by breaking waves on breakwaters is 


the most important thing to be known for the design of breakwater, the problem has been 


~tudied by a number of harbor engineers and investigators during several decades in nature 
as well as in model in many parts of the world. But very few informations have been avail 
able in regard to both the theory and the practical data obtained by measurements in seas, 
because of the difficulties in the experiments and actual measurements due to the high 
intensity and instantaneousness of the shock pressures. Therefore most of structures 
along a shore and in a sea has hitherto been designed by rule of thumb methods established 
from the experiences of disastrous damages caused in maritime structures at storms. To 
the author's knowledge it seems that there exists no generally accepted formula which is 
available with sufficient confidence to insure that both safe and economical designs are 
obtained, being established on the basis of a sufficient number of practical or experimen-= 
tal data. 

In Japan the Hiroi's formula, p = 1.5 WoH, which was established for the pressures of 
breaking waves on the still water surface by Dr. Isamu Hiroi about forty years ago, has 
been used till now for the design of maritime structures by the use of a thumb method, in 
which the pressures of the intensity p = 1.5 w H were assumed to be distributed uniformly 
on the vertical walls of breakwaters from the top to the bottom for the prediction of the 


total wave pressure on the unit length of the vertical walls. 


| 


In other countries the first information was probably given in 1937 by Mr. Larras, 
as far as laboratory data were concerned on the time variation of the shock-pressures of 
breaking waves. ‘r. Bagnold reported in 1939 on the shock-pressures exerted on a vertical | 
wall by a single wave driven forward along a level floor. The actual measurements of ) 


shock pressures were carried out in 1935 and 1937 on the breakwater in the harbor of Dieppe 


But this problem was too difficult to be solved by only these few field and laboratory data 
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The author has been studying on this sibject since 1959 and the experiments are still 
Ww under way. Two thousands or more tests have been performed of the shock pressures 


lerted by breaking waves on the vertical walls of composite-type breakwaters by the use 


two or three wave-pressure-gauges, and the behaviours from instant to instant of the 


‘eaking waves plunging against the breakwaters, the relationships between the wave pres- 
res and the configtrations of the breaking waves, the motion of the water particles in 

e breaking waves, etc. were measured by motion of the 16 millimeter high speed movie taken 
3000 frames per second in synchronism with the pressure gauges. On the other hand, 

ye mechanics of the shock pressures was also studied theoretically. From the results of 

ie theoretical considerations and the experimental data, the new formulas of the maximum 


iock pressure as well as the vertical distribution of maximum simultaneous shock pressures 


sre established for the vertical walls of composite-type breakwaters. 


) Proc. of ASCE, Separate Paper 1063, WW 4, Sept. 1956. 
)} Proc. of ASCE, Separate Paper 1876, WW 5, Dec. 1958. 


) to be published recently in the Proc. of ASCE. 
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| AN EXPERIMENTAL STUDY 
| ON THE SUBMERGED BREAKWATER 


| Masashi Hom-ma, Dr. Engineering 
Professor of Civil Engineering, University of Tokyo 


Toshitsugu Sakou 
Postgraduate Student, University of Tokyo 


Introduction 
A submerged breakwater is a kind of offshore structures, with its top at or below sea 
sr level, used for the creation of a protected area of waters or for the protection of 


stal structures and beaches from damage and erosion caused by wave action. 


It is considered that the submerged breakwater possesses the following two principal 
stions: 
(1) to attenuate waves by causing prematured breaking and partial reflection, 


(2) to bar a seaward movement of the bed materials in the surf zone. 


Up to present there have been a great deal of contributions to the study of these 
stions, theoretically and experimentally, from various sources and from various points 
jiew. Most of them deal with the effects of the height, shape, width, location, etc. of 
breakwater on the transfcrmation of waves passing it. Pre-breaking condition of waves 
sdopted for this treatment. On the other hand, the effect on the sand movement and de- 
nation of beach due to the existence of the breakwater has not sufficiently been argued 

In applying the submerged breakwater to the actual coast, as a means to prevent beach 
sion, however, this latter function of the breakwater might be of more significance and, 
she same time, post-breaking condition of waves should be considered as it is oftener 
ted in the surf zone. Indeed it is our desire to become able to predict the possible 
ige in beach profiles after the construction of the wall and to estimate the amount of 
ir around the wall, for design, construction, and maintenance of the submerged breakwater. 

Therefore, the aim of this study was to know the basic facts in connection with the 
dirensional deformation of beach in an experimental wave flume with breakwaters in place 
especially the depth of scour around the wall placed in the vicinity of the breaking 


It. 


Equipment and procedures 
The wave flume, 17 m long, 0.7 m wide, 0.6 m deep was used, water depth kept at 0.35 m, 
the beach slope 1/15 was buiJt ut cne end of the flume. Fig. 1 shows the mechanical 


ysis of the beach sand used, which was chosen tc reproduce actual conditions and to 
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sec 

1.80 
aig 1.40 220 
11S 170 


Table l 


simplify the analysis by limiting sand size 


variation to a narrow range. 


Table 1 shows the characteristics of 


the model waves employed. A model break- 002 0.1 O02 06 
water was made up of a concrete block stanc : 

and movable wooden blocks which completely Eee 

intercept the flow of water in sand layers b’ neath the breakwater. 

The model breakwater was placed during each run at breaking point, inshore 1/4 of the 
distance from the breaking point to the shoreline, and offshore half of the breaking wave 
length. The ratio of the wall height to the water depth of plecing, h/d, was varied 0.5, 
0.75 and 1.0 for each run and wall location. 


Beach profiles were recorded, for each run, at 15 min., 30 min., 60 min. and 100 min. 


from the start. 


3. Experimental results 

Attacked by waves, initial beach slopes began to ceform in a way corresponding to 
individual wave characteristics, wall location and wall height for each run. Ina fe. 
minutes from the start deformation became visible near the wall, then propagating gradually 
seaward and shoreward. Sand ripples were found throughout the entire bed except for the 
nearest vicinity of the wall. Their size, shape or the region they dominated varied for 
each run and seemed to depend on wall conditions and the curation of the wave attack. 

For small value of h/d, trough-like deformation shoreward of the wall was produced but 
its magnitude was much smaller than troughs formed in case there was no wall in place. It 
is easily imagined that the breakwater acts as an artificial longshore bar in that case. 
For larger values of h/d, there existed the region in the shoreward bed of the wall where 
no visible change or irregularity of bed profiles was observed. 

Figs. 2 and 3 show the change of bed profiles with the lapse of time for wall location 
at breaking point. These figures show that scouring action of waves came to be weak ina 
relatively short time, and then the change of bed profile spreaded gradually seaward. 


Fig. 4 shows the ratio of maxinum scour depth of every 5 minutes to that of 100 mir. 
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from the start for that run. In that figure 
types I, iJ and III represent breaking point, 


inshore end offshore locations of the wall, 


respectively. The depth of wallside scour sea- 


ward of the wall seemed to reach a constant 


value in a comparatively short time. 


Figs. 5,6,7,¢,9 ana 10 show the effects 


of h/d on profile deformation for each set of 


wave and wall location after 100 min. from 


start. Taking Fig. 5, for example, the bed 


immediately seaward of the wall was scuircely 
160 200 240 280 

scoured for h/d=0.50 and 0.75, while for h/d = 
1.00 the depth of scour reached 70 % of the 
water depth of the wall location, and the bar 


formed further seaward grew larger. 


Fig. 6 represents the results with the 


wall in the surf zone, that is, shoreward of 


the supposed trough position without any wall 


in place as is shown in Fig. 11. Comparison il- 


lustrates the decrease of “agnituae of bed de=- 


formation due to the existence of the wall. 
h/d = 0.50 


aie Scour for h/d=1.00 was remarkably deep, but the 


bottom topography is different from that for 
= 10 ie 200 h/d=0.5. Breakwaters placed on a sloping beach 
caused sea level rise shoreward of the wal] as 


ha 2060 they hinder the return flow from the shore. 
=075 


= 1.00 


Flow over the wall due to this rise played an 


important role in scouring shoreward bed, and 


the types of vortice thus produced were dif- 


ferent as is illustrated in Fig. lz. This is 


the cause of the difference of bed scouring pat- 


terns produced by the composite action of waves 


and overflowing waters. 


Fig. 7 is the result with walls offshore: of 


the breaking point. It is noticed that little 


deformation took place around the wall for h/d= 


0.5. That is attributed to the fact that waves 
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sroke spillingly over the wall, exerting little force on the bed in the vicinity. 
Figs. 8, 9 and 10 show similar results corresponding to each of the examples 
pentioned above, depending on each wall location. 
Figs. 13 and 14 are the graphical presentation of the ratios of the scour depth imme- 
ijately seaward of the wall to ‘the water depth of the wall position and the maximum scour 
jepth throughout the entire bed to the water depth of the wall location, respectively. 
From these figures following facts were summarized; 
(1) The ratio is maximum for wall position at breaking point and minimun 
for offshore wall position. 
(2) The ratio increased as h/d increased for breaking point and surf zone 
wall location. It decreased, however, for offshore wall location 


with h/d=1.00. 


(3) The difference of wave characteristics, as far as Wave I and II concerned, 
did not affect the result. 

Summary 2 above can be considered to be due to the supposed fact that the faculty of 
causing scour by overflowing water diminished as the height of the wall increased. 

Fig. 14 employs s(max), the maximum scour depth throughout the entire bed, instead of 
$ in Fig. 13. Difference between s(max)/d and s/d became smaller as h/di increased which 
varifies the explanation on the mechanism of scouring mentioned above. 

Bed immediately shoreward of the wall in Fig. 5 was accreted. Its magnitude increased 
for large value of h/d. Similar statement can be made for Fig. 6. Fig. 7, on the other 
hand, shows that there was no change in the bed vrofile in the shoreward vicinity of the 
wall. As far as all the experiments conducted here concerned, the bed in the shoreward 
vicinity of the wall deformed little and wa: rather accreted. 


ACCRETION HEIGHT 
ON SEA SIDE-FACE 


Bae 
S max * MAXIMUM SCOURING DEPTH AT SEA-SIDE BED 


S SCOURING DEPTH ON SEA-SIDE FACE 
(= MAXIMUM SCOURING DEPTH AT SEA-SID# LED) 


\raatyenes ede 
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Fig. 14 
Fig. 13 = 


Trougn-like nollows produced shore- 


ward of the wall were much smaller than 


these with no walls in place. The sectional 
area of the trough was one third at the 


largest of that with no wall in place. 


It is natural that this magnitude 


depends on the location and the height of 


the wall. Troughs were hardly perceivable 


wnere walls were placed at breaking 
point and large h/d. Offshore wall loca- 


tion did not cause depression of the bed 


near the wall. It moved to the vicinity of 


the shoreline as secondary for breaking of 


waves took place there. 


@ 
234,234 234 
0.50/ 0,75 | 1.00 | 0.50 0.75 


SHOREWARDB.P.| SEAWARDB.P initial position was observed for any of 


Regression of the shoreline from its 


our experiments. This may possibly be due 


Fig. 15 
8 to the fact that the initial slope was much 


steeper than that of ordinary existing beach. 
Fig. 15 shows the ratio of regression for each run to that of no wall in place. Gen- 
erally speaking, submerged breakwaters contributed to lessening the recession of the shore- 


line, but the individual values scattered fairly wide. This may partly be due to sana 


oy OL 
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ples or some irregularity of profiles existing around the shoreline. 


One thing to pay attention to is the rise of water level near the shoreline due to the 


stence of the wall referred in the preceding page, for it makes sea cover more land, thus 
arging the area eroded by waves. Further study is required in this connection. The bot- 
1 Slope of actual beach, moreover, is much gentler and often there is a current along the 
ist which may make the application of the experimental results to actual conditions more 
ficult. Besides, there are other factors such as the concentration of current into the 
s between or the ends of breakwaters which will be serious to planning and maintenance 


yblems. 


Conclusion 

(1) In examining the functions of submerged breakwaters, we must take into account 
their entire effect on beach processes. 

(2) Progress of scour with passage of time depends on the location of the wall and 
the final profile is reached in a relatively short time. 

(3) Seour depth vaies widely according as h/d and wall location vary. The deepest 
relative scour is found for wall location at breaking point and the smallest 
one for offshore wall location. 

(4) Types of correspondence between h/d and the depth of scour for shoreward loca- 
tion differ from those for seaward and breaking point disposition. 

(5) Bed in the shoreward vicinity of the wall was never eroded, but rather accreted. ~ 
This accretion took place rather gradually and indicated the breakwater's 
function as a sand trap against the seaward movement. 

(6) Trough with walls in place were much smaller than those with no wall in place. 

(7) Submerged breakwater evidently decreased the regression of the shoreline but 


we should be careful in its interpretation. 


Studies on beach deformation have much difficulties due to variation and complexity 
factors affecting the phenomena. Construction of an artificial structure will further 
plicate the matter. It is considered that to further the study purposes the post- 
aking characteristics of waves, effects of the bed material properites, a seaward low 
ry the wall and change of orbital motion due to the existence of the wall be evaluated 


completely. In addition, shape effect of the wall must be examined. 
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THE CONSTRUCTION OF A LIGHTED BEACON 
ON A REEF IN THE OPEN SEA 


Yoshio Fujino, Dr. Eng. Director, First Harbor Construction Bureau 
Ministry of Transportation 


ntroduction 

A lighted beacon is one type of aids to navigation erected on isolated shoals or 
s in the sea, Today, there are 125 lighted beacons along the coasts of Japan, but 
are located in port and harbor areas, inland seas, or on elevated foundations such 
mall islands, and it has been considered extremely difficult to construct a lighted 
on on a reef in the open sea, However, recently, we have succeeded in constructing 
acon on a reef off the shore of Omaezaki in Shizuoka Prefecture, 

The reef, known as "Gozeniwa", is located in the Pacific Ocean three kilometers 

the shore, measuring 150 m from the east to the west end, and 1,100 m from the north 
he south end, Though in several sections, the reef stands about one meter above the 
r level, due to the high wind waves and swells in these waters, it is seldom found 
@ exposed above water. 

Past records revealed that very few days of calm weather could be expected to carry 
onstruction works on the site. Moreover, the completed beacon would be exposed to 
vigorous force of ocean waves. 

These circumstances called for special consideration in determining the type of 
eture and method of construction. 

This paper is a brief account of our study of winds and waves in connection with 
construction of the lighted beacon and an outline of the design and construction of 
beacon. 

Study of Offshore Waves Generated by Typhoons 

In order to determine the maximum height and period of waves in the Gozeniwa waters, 
erve as a basis for the design and constructicn of the lighted beacon, we studied 
wind statistics of cutstanding typhoons recorded since the establishment of the 
zaki Weather Bureau. 

Winds of maximun velocity were found among south winds when the center of a typhoon 
ed the west of Omaezaki in the north or northeast direction. 

We studied the records of 583 typhoons which passed this area during the past 67 
s up to 1954, we obtained the location of the center of each typhoon when it was at 
shortest distance from Omaezaki and plotted these positions in quadrants divided by 
lines perpendicular to each other from east to west and north to south with Omaezaki. 
he center of the axis. We numbered the quadrants clockwise beginning with the 


heast quadrant as the first quadrant. We obtained the number of typhoons in each 


THE CONSTRUCTION OF A LIGHTED BEACON ON A REEF IN THE OPEN SEA 
quadrant and classified them according to the distance of the center of the typhoon from 
Omaezaki. Table I gives the frequency of occurrence of these typhoons in the respective 
quadrant and their distance from Omaezaki. 


Table I. Frequency of Occurrence of Typhoons (%) 


reop oneal (ait od 
from Omaezaki (Km 100 200 400 600 800 
Third Quadrant L29 blinded chases 
Fourth Quadrant i lo) alates) 7.4 awe apt: 8 8 BOeS 


Typhoons which are most disastrous to this area are found among those whose centers 


approach Omaezaki in the fourth quadrant. Therefore, we studied 18 outstanding tyhpoons 
which approached Omaezaki in the fourth quadrant since 1891 to obtain their maximum 
velocity and atmospheric pressure. In ten cases, the velocity exceeded that of typhoon 
Tess of 1953, and in five cases, the velocity and the center pressure were comparable to 
those of typhoon Tess. Typhoons of the scale of Tess approached this area five times in 
67 years, which corresponds to 1% of the total number of typhoons and the rate of 
occurrence is once in 14 years. 

Observation of waves had been carried on by the Omaezaki Weather Bureau on the shore 
of Omaezaki from Aug. 1951 to Sept. 1953, and we computed the wave height of offshore 
waves from the observation data. 

Waves of maximum observation values were those generated by typhoon Tess, with a 
period of l2 second:. Outstanding typhoons recorded during this period include typhoons 
Marge, Ruth, Dinah, Jean, and Nina. 

We computed the height and period of waves generated by Tess and the other five 
typhoons according to the Sverdup-Munk method, the results of which are presented in 
Table 2, 


Table 2. Waves of Maximum Dimensicns in the Omaezaki Waters 


Date 


Sept. 1934Joct. 1949]0ct. 1951]Sept. 1953] Sept. 1954 | 


2. Name of Typhoon 


Lorna 
Maximum Center pré@8ure 


Minimum distance of 
center from Omaezaki 


Direction of maximum 
wind at Omaezaki 


6. Maximum velccity of 
wind at Omaezaki(m/sec) 


ee 


a ee 
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Jinds which generated 
aximum waves 

7. Velocity m/sec 
8. Fetch (knot) 

9. Duration (hour) 


lO. Distance of decay 
(knot) 


Vaves of maximum dimensions 
ll. Wave height of 
significant waves (m) 


l2. Period of significant 
waves (sec) 


13. Travel time 


When the Bretschneider and Piersen-Neumann-James methods are applied to typhoon Tess, 
c-esults are as given in Table 3. 


> 3. Waves of Maximum Dimensions generated by Typhoon Tess 


Basic Theory Sverdup-Munk [Bretschneider) Piersen-Neumann James 


Significant Waves 
9.0 9.2 
8.0 © 


There is a considerable difference among the periods obtained by the three methods, 


Wave height (m) 


Wave period (s) 


the Bretsneider method seems to be reasonable and the calculated value is in close 
sment with the period 12 seconds obtained from actual observation. 

If we assume the maximum period of waves generated by typhoons to be 12 seconds, we 
sbtain the conditions under which these waves break by applying Iversen's chart, the 
't of which is given in Table 4. 

2 4, Depth of water and wave height of breakers 


Wave height of 
offshore waves Ho (m) 


Wave period of 
offshore wave T (sec) 


Wave height of breaker 
Ho (m) 


Depth of water hb (m) 


The waves break in the vicinity of Gozeniwa where the water depth is over 10 m and 
ave height reaches 80 - 8.7 m. However, as the water depth within the Gozeniwa 

area is estimated to be 9 m, offshore waves do not directly reach this area. Waves 
is area would retain the maximum period of offshore waves, but the wave height 


pe limited by the water depth. Table 5 gives the maximum dimensions of waves 


= Sis 


THE CONSTRUCTION OF A LIGHTED BEACON ON A REEF IN THE OPEN SEA 
breaking in waters with depths of 9 m obtained from Iversen's chart. These waves may 
be considered to be of the maximum height within the reef area. 

Table 5. Waves of Maximum Dimensions in Water Depths of 9 m 


Water depth hb (m) 
Period T (sec) 


height of offshore 
waves Ho (m) 


height of breakers 
Hb (m) 


III. Investigation of Waves in the Reef Area 

In view of the difficulty in analyzing the condition of waves within the Gozeniwa 
reef area, we made a qualitative analysis of the characteristics of waves generated by 
typhoons such as Tess which attacked Gozeniwa from the southern areas at 1 P.M. April 


18, 1956 by means of aerial photographs. (See Fig. 1) 


Fig. 1 Aerial Photograph of Waves in the Gozeniwa Waters 


The following facts were revealed by the photographs. 
(a) Offshore waves from the southern seas interfere with each other before they 
reach the Gozeniwa reef and no longer possess the form of offshore waves at the 


northern end of the reef. However, they maintain the wave period of offshore 


waves. 
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(b) From the reflection of the sun rays, it may be assumed that the wave height is 
| highest at the southern end and decreases as it travels northward. 

(c) Around the northern end of the Gozeniwa reef, there is an interference of waves 
which may produce choppings. 


(a) Therefore it would be advisable to select a locality in the northern part of tne 


reef where the wave force has been reduced as the construction site of the 
lighted beacon. 
- Wave Height and the Design of ‘the Lighted Beacon 

The lighted beacon was decided to be constructed on the northern end of the reef in 
sw of the above factors. The maximum height of waves in this locality depends on the 
pth of water. However, in this case, we must consider the tidal range and extraordi- 
ry rise of water as well as the average water depth.of the area. 

The water depth at this locality is -6 m on the chart and the tidal range is 2.0 m. 
times of typhoons, assuming that the decrease in atmospheric pressure would be 100 
S, and the extraordinary rise of water would be 10 m at the most, the maximum water 
oth at the site of construction would be 9.0 m. 

If we determine the wave period of maximum waves to be 12 seconds according to 
ble 3, the maximum wave height in this locality would be Hb = 7.8 m. This figure was 
opted in the design of the substructure of the lighted beacon with due consideration 
r the interference of waves and choppings generated by typhoons. 

Construction Works and Meteorological Conditions 

The number of calm days which would be available for construction works on the sea 
re studied as a basis of determining the design and method of construction. 

The standard for determining whether it is possible or not to carry on the operation 
works on the sea differs greatly according to the type and nature of work and the 
sign of the structure, and Glenn of America has carried on an investigation in this 
91d. According to his studies, in regard to works on lighted beacons, generally, wave 
ights of 2 ft or 0.6 m mark the limiting point of conducting operations on the sea. 
we compute the wind velocity in this case from a formula for the relationship of wind 
locity and wave height, V = 3.3 m/sec, belonging to Grade 2 in the wind force order. 

We calculated the number of calm days to be expected for construction works on 
3 sea by the following methods. 

1. The number of calm days available predicted form the duration of weak winds. 

Wind velocity and wave height are decisive factors in determining whether the 
conditions are suitable for works on the sea or not. However, as data on the 
wave height at Gozeniwa were not available, we calculated the number of calm days 
from the duration of weak winds and obtained results as shown in Table 6. Most 


favourable conditions were found in July and August, but it is extremely 
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difficult to expect calm weather to continue over 48 hours. 


Table 6, Duration of Winds below Wind Force Order 2 (3.3 m/sec) 
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2. The number of calm days available determined from the analysis of wave height and 
wind velocity. 

In order to determine the number of calm days available from the point of wave 
height as well as wind velocity, we analyzed the statistical data from observa- : 
tions carried on during the past 4 years, at 24 lighthouses throughout the country, 
We picked out the calm days when the wind waves, swells and wind velocity 
belonged to Grade 1 or 2 in each water zone and season. On the Pacific coast, 
the number cf calm days when both wave height and wind velocity were in favorable 
condition was 70% of the number of calm days determined by wind velocity alone. 

Hence, with reference to Table 6, we may forecast the number of calm days 
available for construction works in the Gozeniwa waters. 

On the other hand, if it is assumed that a wave height of 0.5 m recorded in 
the wave observation in the Omaezaki harbor is the maximum limit for the 
operation of works on sea, the maximum wind velocity in this case would be 7.9 
m/sec, slightly stronger than desirable. Calm days with wave heights and wind 
velocity below these values are given in Table 7. In July, there are the largest 
number of calm days, the maximum being 10 days, and May, June, September and 
April follow in order with the number of calm days. In August the number of calm 
days varies greatly each year due to the attack of typhoons in this month. 

Table 7. Number of Days of Calm Seas (Max. wind velocity 7.9 m/s) 


(Max. wave height 0.5 m ) 
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3. Forecasting the number of calm days from the weather chart. 


ia Se : 
3.6 


Selecting cases when favorable wind velocity and wave height continued for 
more than 48 hours, we studied the pattern of the weather chart in each case and 
obtained the following conélusions. 

(a) Construction works on the sea should be carried on during the period from 
April to September, of which July will be the most favorable month. 
(bh) As a rule, favorable weather may be expected when variable high atmospheric 
pressure or the north Pacific high atmospheric pressure overlies the area. 
(ec) The cold front extending from the south to north, and high atmospheric 
pressure traveling at a high velocity are apt to bring unfavorable weather. 
. The Substructure of the Gozeniwa Lighted Beacon 
The substructure of the lighted beacon to be constructed in the open sea was bound 
' various restrictions from the standpoint of designing and method of construction. 

(a) The structure must be of sufficient strength to withstand the vigorous force 
of ocean waves. 

(b) The principal part of the construction works must be carried on at the work 
base so that works on the sea will be limited and simple. 

(ce) Works on the spot must be completed in a short period of time. 

(ad) Each unit of work on the spot must be completed within one round of continuous 
calm weather, and the completed portion must be able to withstand the attack 
of minor storms until the next round of continuous calm weather when works 
will be undertaken again, gradually increasing the stability of the structure 
with each unit of work. 

A braced cylinder type tripod structure as shown in Fig. 2 was adopted as 


the substructure to meet the above conditions. 
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Fig. 2 The Gozeniwa Lighted Beacon 
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The body, with a platform, was assembled at the base, launched in water and towed 
e site of construction where it was fixed on the reef 6 m under the water level, 
rovided to withstand swells of 1 meter. 
Then I type steel piles were driven into the reef through the cylinders to stabilize 
tructure to withstand swells of 2 - 3 meters. 
The cylinders were then filled with concrete, and prepact concrete was cast around 
hree cylinders at the base to form a triangular body, one side of which measures 
m to a height of 1.8 m, firmly adjoining the structure to the bed rock. 
Theories of Wave Pressure exerted on Lighted Beacons 
Theories of wave pressure exerted on submarine structures may be divided into two 
categories, wave pressure on vertical walls, and wave pressure on cylindrical 
S, Wave pressure exerted on lighted beacons may be classified in the second category. 
Theory of wave pressure on cylindrical bodies 
Recently, with the development of new theories, various experiments have been 
conducted on computing the wave force acting on cylindrical bodies in water. 
Generally, friction drag F, and inertia drag Fp act on dg a portion of a 
cylindrical body in fliud. 


When these two are represented by F. 2 


£ 
= = 2 D du rcterslee cig cise ntee ase) 
Cp : drag coefficient u : lateral velocity of water particles 
Cy : inertia coefficient Ff: specific weight of water 


D : diameter of cylindrical 
body 
The coefficients Cp, Cy vary according to the motion of water on the cylinder, 
and may be considered to be constant for each part of the cylindrical surface. 
The maximum moment does not occur at the wave crest but at a point with an 


angle displacement from the crest. 


= aint / RM D (aac elt ELE Seine cena) 
8 Cp BH Kop=. (0 4:2Z)_K3} 

Ky, Ko, K35 K, are functions of Z only. 

When the water depth d is relatively small compared to the wave length L, or 
when the diameter of the cylinder D is small compared to the wave height H, the 
value of B approaches zero. On the other hand, when d/L and D/H are large values, 
the value of B approaches 90°. The values of Cp, Cy are determined by the wave 
curve and moment-time curve. 

Experiments of J. R. Morison 


Morison erected a cylindrical structure in water and conducted observational 
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experiments on offshore waves, swells starting to break, breakers, and breakers 


with some foam, and obtained the values of the drag coefficient Cp (See Fig. 3) 
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Fige 3 Variation of Cp in relation to H/d 
The experiments were conducted under the following conditions: 
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Therefore D/H and d/L were both small values and may be assumed to be almost 
zero. 

Hence Cy = 0, 

and B was not actually obtained in the experiment. 

According to the diagram, the value of Cp is greater in the case of breakers 
compared to swells, indicating a value close to 1.0. Strictly speaking, this 
theory is not applicable to breakers, but actually, the results are in close 
agreement with the general conditions when breakers act on lighted beacons, and 
it may be considered to be an outstanding experiment in this field. 

Resistance force of cylindrical bodies in steady flow 
Generally, the pressure exerted by waves on cylindrical bodies in steady flow 


may be represented by the first term in equation (1) from which the drag coeffi- 
cient C, may be obtained. Cy is determined by the Reynold number (Re) only, and 
e ste f 


it shows a sharp decrease and reaches a transient state when, 


Re = 3-5 x 10? 
but enters a steady state after it passes this point. 
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Fig. 4 Drag Coefficient of Cylindrical -Bodies in Steady Flow 
When we compute the value of the Reynold number in the case of J. Re Morison's 
experiment, the value obtained is within the transient state, and a wide scatter 


may be observed in the value of Cp: 
In the case of lighted beacons, we apply the formula Re, = 10° 107, and the 


value of C, may be assumed to be in a steady state. 
Lateral force exerted by waves on conerete blocks on the sea bed. 

T. Hamada has conducted experiments on the pressure exerted on concrete blocks 
on the sea bed in shallow water and breaking points. 

In the case of swells, the inertia coefficient is of a considerably large value 
and in the case of breakers, it is 20% - 40% of the drag coefficient. When x 
is the length of one side of the block, 

so fol = OHS; — IO) 
However, in the case of lighted beacons, it is generally, 
D/H = 0.1 - 0.2 

If a similar experiment were conducted on a lighted beacon, the inertia 
coefficient would most likely be of a negligible value. 
Computation of the wave force on pier piles 

Munk has computed the wave force on a cylindrical pier pile off the shore of 
Scripps Institute. He has computed the wave pressure exerted by breakers by the 
solitary wave theory with Cp = 0.33. 

However this computation involves some questionable points. 

(a) As the behavior of breakers differ from steady flow, even when Re=~ 5 x 107, 
the drag coefficient should be greater than Cp = 0.33. 

(b) Though Cp = 0.33 may be applicable to the drag coefficient at the crest of 
the waves, as the water particle velocity distribution varies along the 
vertical surface ‘here is some question in giving a constant value for Cp. 


(c) There should be some consideration for the sharp increase in wave height 
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when breakers or waves impinge on the front of the structure. 
6. Theory of R,Curtis Cooke 

Cooke applied the results of Iversen's investigations and proved that the values 
of Cy and Cy differ for each wave phase. He proposed a common drag coefficient 
C and proved that C may be obtained from the Iversen Modulus and the Reynolds 
number. 

However, when we compute the drag coefficient of breakers by the solitary 
wave theory, the value of the Iversen modulus will be zero, as the water particle 
velocity reaches a maximum, and the acceleration will be zero. Hencé the drag 
coefficient will be determined by the Reynolds number alone. 

7. Summary of theories of wave pressure on lighted beacons 

From the study of various theories and experiments on wave pressure exerted on 
cylindrical bodies, we have obtained a satisfactory approximation to the wave 
pressure exerted by breakers on lighted beacons. 

The inertia coefficient may be disregarded. On the other hand, it has been 
proved through experiments that the vertical distribution of wave velocity at the 
breaking point closely resembles that at the limit of breaking according to the 
solitary wave theory. R. L. Wiegel has also reported similar results from his 
experiments, 

Therefore, the pressure exerted by waves on lighted beacons may be said to be 
determined by the resistance force of a cylindrical body exposed to the pressure 
of solitary waves at the limit of breaking. 

When the pressure exerted by waves on a lighted beacon is Fo, 

Fo = $P Cy fdu2 az 
With reference to the experiments of J. R. Morison, allowing a sufficient 
margin for safety, we assumed that, 
Cy = 1.0 
in our computation. 
VIII. Conclusion 

I have presented a brief account of the design and construction of a lighted beacon 
on an isolated reef in the Pacific Ocean. 

However, in regard to the vertical wave pressure exerted by choppings on submarine 
structures, I have merely presented problems yet to be solved, and there is also much to 
be explored on theories on wave pressure on cylindrical bodies. 

Nevertheless, I hope that our technical approach to the problem of planning and 
constructing a structure in the open sea, stepping into a hitherto undeveloped field, 


will serve to promote the development of engineering in this field in the future. 
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